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INTERIM  FEASIBILITY  ASSESSMENT 
METHOD  FOR  SOLAR  HEATING  AND 
COOLING  OF  ARMY  BUILDINGS 

1  INTRODUCTION 

Background 

Like  all  energy  users,  the  Army  is  faced  with  a  rap¬ 
idly  increasing  energy  bill  and,  in  some  locations,  a 
shortage  or  curtailment  of  its  energy  supply .  Therefore, 
the  Army  must  look  for  new  natural  energy  sources, 
especially  those  which  are  abundant  and  are  not 
subject  to  inflation.  The  most  ideal  so,.rce  is  solar 
energy,  because  of  its  abundance,  widespread  distribu¬ 
tion,  and  absence  of  recurring  fuel  cost.  Therefore, 
solar  energy  systems  offer  tremendous  potential  as 
alternate  energy  sources  fur  heating  and  cooling  Army 
installations,  if  technically  feasible  utilization  systems 
can  be  developed  which  are  economically  competitive 
with  conventional  fuels. 

The  technical  feasibility  of  heating  and  cooling 
buildings  using  flat-plate  solar  collectors  has  been 
established  both  in  theory  and  practice.  Although 
future  demonstrations,  including  the  Army’s  plans  for 
a  solar  demonstration  at  Port  Hood, TX,  will  indicate 
improvements  in  component  design  and  manufacturing 
and  in  system  design  methods,  investigators  can  ap- 

iiuip.!  I  1 IV-U I  il  mIIO  C«~II1II^  1  fV  1 1 1 1  iTfli 

confidence  th':t  a  practical,  reliable  system  can  be 
constructed.  Although  the  design  phase  may  be  some¬ 
what  more  complex,  the  construction  phase  requires 
little  nKirc  skill  than  is  required  to  install  conventional 
‘leafing  tmd  emsiing  systems.  Tlie  decwtmi  whethet  tn 
apply  solar  energy  fur  heating  and  cooling  should 
therefore  generally  be  based  on  life  cycle  cost  compar¬ 
isons,  takmg  into  consideration  the  escalating  price 
and  decreasing  availability  of  conventional  fuels. 

Purpose 

The  purpose  of  this  report  is  to  provide  a  method  to 
p«t4<  m  pnfimlnuy  eteonssmte  feeslMfty  autsssmetM 
of  candidate  solar  heating  and  cooling  systems  when 
applied  to  specific  buildings  at  specific  sites. 

Approach 

The  economic  assessment  of  solar  energy  systems 
mtohm  cpinpaThig  eonrwrtfanal  Wt-cytk  ttsel  costs 
with  the  additional  first  cost  of  installing  a  solar  energy 


system  and  the  resulting  reduced  operating  costs.  The 
principal  problem  of  tliis  analysis  is  determining  the 
amount  of  energy  that  the  solar  energy  system  can 
economically  supply  to  meet  the  heating  and  cooling 
requirements  of  specific  buildings  and  sites.  Unfortun¬ 
ately,  simple .  design  calculations  generally  do  not 
provide  a  basis  for  determining  either  the  annual  build¬ 
ing  energy  requirements  or  the  fraction  of  those 
requirements  that  can  be  met  economically  wi'h  solar 
energy.  Computer  simulation  techniques  r.iust  usually 
be  used  for  these  design  calculations. 

To  develop  the  leasihility  assessment  techniques, 
computer  simulation  was  used  to  perlorm  detailed 
analysts  ol  mmieroiis  candidate  solar  heating  and  cool¬ 
ing  systems  lor  test  case  buildings  at  vaiious  sites. 
Based  on  analysis  results,  a  general-purpose,  dimension¬ 
less  curve  was  developed  which  can  be  directly  applied 
to  other  buildings  and  sites  to  estimate  expected 
performances  of  candidate  solar  heating  and  cooling 
systems.  Using  this  estimated  iierformance,  a  simplified 
method  was  developed  to  initially  assess  the  potential 
economic  viability  of  the  solar  project  without  expen¬ 
ding  considerable  funds  and  effort  for  detailed  analysis. 

Scop* 

Only  heating  and  cooling  applications  using  solar 
energy  were  considered,  so  assessment  results  do  not 
apply  when  only  solar  heating  is  being  considered. 
ConMne^  «oIm  beM^ig  and  «iMlf(n|  mu 
consideration  primarily  because  the  annual  energy 
requirement  fur  cooling  heavily  dominated  that  of 
heating  for  the  buildings  and  sites  studied.  This  domi¬ 
nance  can  be  expected  as  long  as  current  comfort 
«ei>ISng  deal*!)  ttenda  -Me  appHed. 
for  solar  heating  only  are  presently  being  developed, 
and  these  will  be  the  subject  of  a  separate  report. 

Chapter  2  of  (his  report  brielly  ('escribes  solar  heat¬ 
ing  and  cooling  systems  and  their  components.  Chapter 
3  presents  the  simulation  studies  results  and  the  curves 
and  equations  necessary  to  estimate  expected  solar 
heating  and  cooling  system  performance  for  buildings 
atwl  a  a 

method  for  determining  the  lowest  life-cycle  cost  solar 
heating  and  cooling  system  and  for  determining  its 
economic  viabiliiy.  l^apief  presents  ihe  report’s 
conclusions  and  recommendations. 

Appendix  A  provides  an  example  application  of 
rtie  twiiomic  fetsIbiVrt^  irssesstneiYt  descrditd 

in  Chapter  4.  Appendix  B  describes  the  National 
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Bureau  ol  Standards  Load  Determining  Program 
(f^BSLD)  application  to  the  barracks  and  battalion 
headquaiters  and  classroom  buildings.  Appendix  C'  is 
a  detailed  description  of  the  solar  heating  and  crnrling 
sinuilatioM  model. 

2  COMPONENTS  OF  A  SOLAR 
INSTALLATION 

Gtnaralixad  Syttam 

Figure  I  is  a  schematic  diagram  of  a  basic  system  for 
heating  and  cooling  with  solar  energy.  In  comparison 
to  conventional  systems,  the  only  unique  features  are 
the  collector  array  and  the  thermal  energy  s' orage  tank. 

The  sunlight  falling  un  the  array  warms  a  fluid 
(usually  glycol  and  water),  which  is  pumped  through 
the  solar  collectors.  The  heat  from  this  fluid  is  trans¬ 
ferred  within  a  heat  exchanger  to  a  second  fluid  (usually 
water),  which  is  pumped  to  an  insulated  thermal  stor¬ 
age  tank.  When  there  is  a  heating  or  cooling  demand, 
the  warm  fluid  from  the  storage  tank  is  either  pumped 
directly  to  the  heat  exchanger  in  the  duct  or  to  the 
cooling  unit.  In  most  feasible  systems, the  temperature 
in  the  tank  will  occasionally  drop  below  useful  values; 
it  is,  therefore,  necessary  to  include  in  the  system  an 
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auxiliary  heater  capable  of  supplying  all  or  part  of  the 
heating  or  cooling  demand. 

Solar  Collacton 

The  function  of  a  solar  collector  is  to  capture  as 
much  available  sunlight  as  possible  and  to  convert  this 
energy  to  heat  energy.  Converting  solar  energy  to  heat 
at  low  and  moderate  temperatures  (80°F  to  300°F 
(24.6“C  to  l48.7°Cp  involves  components  which  are 
in  low-volume,  commercial  production.  This  tempera¬ 
ture  range  is  adequate  for  heating  and  cooling  buildings 
and  will  be  considered  in  this  report. 

Two  types  of  solar  energy  fall  on  a  collector:  that 
coming  directly  from  the  sun  (beam  energy)  and  that 
coming  from  other  directions  through  scattering  or 
rcfleclk)n  of  the  sun's  rays  (dilTuse  energy).  On  a  clear 
day,  the  intensity  of  diffuse  energy  is  much  less  than 
that  of  beam  energy,  but  the  former  is  often  intense 
enough  to  be  useful  during  cloudy  or  smoggy  condi¬ 
tions. 

Tracking  and  Concentrating  Collacton 

To  maximize  the  amount  of  beam  energy  collected, 
the  solar  collectors  can  be  continually  moved  so  that 


they  renuin  oriented  toward  the  sun;  i.e.,  tracking  the 
sun.  The  collected  energy’s  density  can  be  further 
increased  by  using  focming,  or  otliei  concentrating 
methods  to  intensify  the  sun's  radiation  on  the  collec¬ 
tors.  Although  this  can  result  in  heat  of  higher  temper- 
atuHii,  iTwant  for  irM'kiug  t’ld  ctirwnnulitng  are 
presently  expensive.  Consequently,  this  report  will 
consider  only  fixed,  flat -plate  collectors. 

Flat -Plata  Collsctort 

Figure  2  is  a  cutaway  sketch  of  a  conventional 
flat-plate  collector.  The  unit  is  mounted  in  a  fixed 
position,  and  the  metallic  (copper,  aluminum,  or  iron) 
absorber  plate  is  coated  to  maximize  solar  energy 
absorption.  Often,  the  coating  has  low  emissivity  in  the 
infrared,  which  minimizes  the  re-radiation  of  energy 
from  the  hot  absorber  plate.  Heat  loss  to  the  environ¬ 
ment  caused  by  conduction  through  the  back  of  the 
collector  is  kept  low  by  heavy  insulation.  The  trans¬ 
parent  cover  plate  reduces  heat  losses  from  convection, 
and  by  having  low  transmittance  to  infrared,  reduces 
radiation  losses.  This  design  enhances  the  “green¬ 
house”  effect  and  maintains  a  plate  temperature  which 
is  much  greater  than  ambient  air  temperature.  Some 
models  use  two  cover  plates,  separated  by  an  air  space, 
to  Bmhet  reduce  cumccHun  luss,  tftoUgFi  thh  Is  at  the 
expense  of  transmittance. 

Useful  heat  is  extracted  from  the  collector  by 
circulating  air  or  a  liquid  over  or  through  the  absorber 
plate.  Under  fav»)rable  conditions,  the  heat  energy 
removed  can  exceed  50  percent  of  the  solar  energy  in¬ 


cident  on  the  panel.  For  a  given  collector, el ficiency  is 
highest  when  the  plate  is  coolest  and  when  the  differ¬ 
ence  between  plate  and  amlitent  lempeTaime  is  mini¬ 
mal. 

Uigit  IhcmWhi  ioMi  PwlwlhMt  sttd  llw  alMMrk.>v 
collector  fluid  circulation  (stagnation)  can  cause  the 
collector  plate  temperature  to  exceed  .150 'F'  (I77"(’). 
Thus,  enritector  materials  ntust  be  selected  to  withstand 
a  wide  temperature  range  and  high  thermal  gradients. 

.Several  manufacturers  market  solar  colicclois  with 
wide-ranging  collector  performances  and  costs.  I  he 
National  Bureau  of  Standards  has  puhlishid  a  manual, 
NBSIR-74-().15,  which  outlines  test  procedures  designed 
to  permit  collector  performance  comparisons.' 

Colltctor  Fluids 

The  most  common  fluids  used  to  transfer  heat  from 
the  collector  are  air  and  water  or  water-glycol  solutions. 
Liquid  systems  are  generally  preferred  to  air  for  com¬ 
bined  heating  and  cooling  systems,  because  less  power 
is  required  for  circulation;  there  are  fewer  problems 
with  transmission  and  insulation  (ducts  vs.  pipes);  there 
is  smaller  storage  mass  and  volume;  and  there  is  no 
commerciaDy  availabte  cuoting  riiethod  width  tan 
directly  use  heated  air. 


'j.  S.  Hill  and  T.  Kusuda,  Methods  of  Testing  for  Ratinn 
Solar  Collectors  Based  on  Thermal  Performance,  NBSlR-74-635 
(Thermal  Kngineering  Systems  Secliun,  Center  for  Building 
Technology,  National  Bureau  of  Standards,  December  1974). 


Figure  2.  Conventional  flat-plate  collector. 


Wutcr  muy  be  circulated  llirouch  collector!!  in 
locations  wliere  f'rce/ing  docs  not  occur,  or  in  colder 
regions  if  the  system  is  drained  and  Hushed  with  inert 
gas  whenever  sulTicicnt  solar  energy  is  not  available; 
however,  the  must  commonly  used  liquid  is  a  water- 
glycol  solution,  bulTercd  to  reduce  corrrrsion.  In  addi¬ 
tion  to  decreasing  the  water's  I'rcc/ing  point,  the  glycol 
increases  its  boiling  point.  This  is  an  important  advan¬ 
tage,  since  when  the  collector  is  operational  its  temper¬ 
ature  can  rise  above  2I2°F  ( l()0°(').  A  relief  valve  and 
an  expansion  tank  must  be  included  in  the  collector 
system  to  insure  its  safe  operation. 

To  reduce  electrolytic  corrosion,  both  within  the 
collectors  and  elsewhere,  piping  ol  diffwrent  metals 
must  he  isolated  by  fittings  or  by  piping  which  is  not 
electrically  conductive.  If  a  water  solution  is  used,  the 
pH  must  be  monitored  and  buffers  and  inhibitors 
added  periodically  to  suppress  corrosion. 

Other  low-viscosity  heat-transfer  liquids,  c.g.,  sili¬ 
cones,  may  deserve  consideration  as  collector  system 
Huids,  because  they  can  attain  higher  boiling  points 
and  possibly  minimize  corrosion  problems. 

Colltctor  Location 

The  most  common  location  for  an  array  of  collectors 
is  the  roof  of  the  building  to  be  heated  or  cooled.  If 
the  roof  slope  does  not  coincide  with  the  desired  col¬ 
lector  tilt  angle,  the  collectors  shotild  be  spaced  to 
insure  that  one  row  does  not  significantly  shadow 
another  at  any  time.  If  several  buildings  must  be  served, 
or  if  the  roof  area  is  too  small,  it  may  prove  feasible 
to  mount  collectors  on  an  unobstructed  land  area 
adjacent  to  the  buildings  being  served.  The  collector 
array  should  generally  be  as  close  as  possible  to  the 
energy-using  facilities  to  minimize  piping  costs  and 
heat  loss,  and  to  reduce  pumping  power  requirements. 

Collector  Araa 

Since  the  cost  of  the  collectors  is  a  substantial 
portion  of  the  total  solar  installation  cost,  the  number 
of  collectors  (or  tiie  total  useful  collector  area)  must 
be  carefully  determined.  The  required  collector  area 
will  be  proportional  to  the  load  being  met  and  to  in¬ 
cident  solar  radiation,  both  of  which  vary  according  to 
site.  Chapter  3  discusses  determination  of  the  collector 
area  based  on  computer-derived  curves. 


Thermal  Storage 

Heat  storage  capacity  is  crucial  in  a  solar  installation 
design.  Since  the  incident  radiation  on  the  collectors 
is  only  enough  to  be  useful  for  a  few  hours  each  day,  it 
is  necessary  to  store  heat  ■'j  thai can  accommodate 
part  or  all  of  the  demanded  load  a.  night  or  during 
cloudy  weather.  Ideally,  it  would  be  desirable  to  store 
heat  during  a  higli  radiation/low  load  period  for  delivery 
during  a  low  radiation/high  load  period,  i.e.,  season  to 
season.  Cost  and  size  presently  make  such  large  storage 
capacities  marginal,  but  storage  capacity  to  handle  one 
week's  load  now  seems  feasible. 

Solar-dcrivcd  heat  may  be  stored  either  as  sensible 
heat  or  as  latent  heal.  Latent  heat  is  used  to  melt  a 
solid,  such  as  (ilauher  salts  or  paraflln.  Although  such 
systems  have  small  volume  and  mass,  high  cost  and  low 
reliability  presently  limit  their  use.  Water  is  the  best 
means  of  storing  sensible  heat,  primarily  because  of  its 
lower  cost.  To  be  useful,  the  water  temperature  must 
exceed  the  temperature  demanded  by  the  load-approx- 
imatcly  90“F  (26.4°C)  for  heating,  and  ISO^F  (81 .4“C) 
for  absorpw:./n  cooling.  The  upper  temperature  limit  is 
the  boiling  point  (212‘’F  [100°C]  at  sea  level).  The 
tank  may  be  made  of  metal,  concrete,  or  fiberglass, 
depending  on  which  is  locaUy  most  economical.  The 
surface-area-to-volume  ratio  can  be  minimized  with  a 
cylindrical  tank  whose  diameter  is  equal  to  its  height; 
for  a  tank  with  plane  sides,  a  cube  is  best;  a  sphere  is 
ideal. 

The  storage  tank  must  be  well-insulated  so  that 
heat  loss  through  its  walls  will  be  very  small  in  compar¬ 
ison  to  the  heat  delivered  to  the  load.  Since  the  surface- 
area-to-volume  ratio  becomes  smaller  with  increasing 
tank  size,  the  required  insulation  thickness  is  less  for 
larger  systems.  Modern  insulation  materials,  such  as 
polyurethane  foam,  provide  minimal  heat  loss  at 
moderate  cost. 

The  NBS  report.  Method  of  Testing  for  Rating 
Thermal  Storage  Devices  Based  on  Thermal  Perfor¬ 
mance,  provides  a  method  for  evaluating  storage  sys¬ 
tems.^ 


*G.  E.  Kelly  and  J.  E.  HWi,  Method  of  Testing  for  Rating 
Thermal  Storage  Devices  Based  on  Thermal  Performance, 
NBSlR-74-634  (Thermal  Engineering  Section,  Center  for 
Building  Technology,  National  Bureau  of  Standards,  May 
1975). 


10 


Htating 

Conventional  means  arc  used  to  heat  a  building  with 
stored  solar  energy.  Hot  water  from  the  storage  tank  is 
pumped  through  a  coil  in  the  air-handling  system, 
where  the  air  is  heated  and  delivered  to  the  building. 
Unlike  a  conventional  system,  however,  the  water 
temperature  varies  considerably  (80"F  to  2I0“F 
|24.6‘’C  to  81 .4°C| );  therefore,  the  fans,  pumps,  and 
coils  must  be  capable  of  providing  good  heat  exchange 
at  the  lower  temperatures.  The  wide  range  of  water 
temperatures  makes  radiators  and  baseboard  convectors 
less  attractive  than  forced,  central  air  systems.  Heat 
stored  in  the  tank  can  also  be  used  to  heat  or  preheat 
water  fur  domestic  use  in  the  building. 

When  the  tank  lemperaluie  drops  below  the 
minimum  useful  temperature,  (usually  about  80“F 
|2‘>.(i"C|),  controls  must  be  provided  to  turn  on  the 
auxiliarv  heater,  which  is  usually  a  conventional  gas. 
oil,  or  electrical  heater.  .Since  this  auxiliary  heat  is 
expensive,  controls  and  piping  should  be  provided  so 
that  the  wal'ii  bypafsaes  the  ihiwagc  lartk  when  (he 
auxiliary  heal  is  in  use. 

Alternate  heating  methods  may  be  used.  A  heal 
pump  can  provide  heating  at  considerably  lower  tank 
temperatures;  however,  since  solar  energy  cannot  be 
used  to  drive  the  heat  pumps  for  cooling,  this  method 
IS  suitable  only  for  solar  heating  applications. 


l  ithium  bromide  and  water  vapor  absorption  cool¬ 
ing  units  which  were  originally  designed  for  gas,  steam, 
or  hut  water  firing,  have  been  developed  to  operate 
at  generator  temperatures  as  low  as  I8()''F  (81.4 ’(') 
with  a  coelTicient  of  performance  (COP)  belter  than 
These  units  are  commercially  available  in  limited 
si/es  and  have  been  successfully  operated  in  solar 
installations.  Manufacturers  of  large-capacity  absorp¬ 
tion  chillers  are  beginning  to  examine  the  potential 
applicability  of  their  equipment  to  solar-powered 
systems.  Absorption  cooling  is  currently  the  recom¬ 
mended  method  for  cooling  with  solar  energy,  because 
it  has  been  successfully  demonst rated,  and  because  it 
is  currently  the  most  inexpensive  solar  cooling  method, 

Raliibility 

The  design  of  a  solar  installation  must  consider  both 
reliability  and  maintenance  cost.  Since  the  solar  collec¬ 
tors  and  the  thermal  storage  tank  are  essentially 
Matk,  mafnleflanue  costs  and  system  lehability  sItuuUi 
be  comparable  with  those  of  conventional  heating  and 
cooling  systems.  An  NBS  manual  prepared  for  HUD, 
Intenm  Performance  Cniena  for  H)tar  Healing  and 
Owling  Systems  and  Dwellings,^  provides  some 
reliability  criteria. 

3  SYSTEM  PERFORMANCE  ESTIMATES 

Simulation  Study  Prooadurt 


Cooling 

Cooling  a  building  by  solar  energy  can  be  accom¬ 
plished  by  several  methods,  including  Rankine  cycle, 
desiccant,  and  vapor  absorption  cooling. 

Solar-powered  Rankine  cycle  turbine  engines  have 
been  used  to  drive  conventional  vapor  compression 
ciKiling  systems'  however,  these  systems  are  experi- 
mfrvtat  and  piuliibilitoty  ovpenkiw,  to  they  ah  not 
now  considered  to  be  a  practical  alternative. 

A  recently  introduced  prototype  commercial 
desiccant  system,  designed  to  be  powered  by  solar 
energy  and  natural  gas,  performs  particularly  well  in 
dry  climates;  however,  its  performance  coefficient  in 
humid  regions  may  be  lower  than  that  of  absorption 
coolers.  Since  this  system  is  not  yet  commercially 
available,  it  will  not  be  considered  in  this  report. 


In  the  past,  selection  or  sizing  of  conventional 
heating  and  cooling  systems  was  accomplished  with 
relatively  straigntforward,  steady-state,  peak-load  cal¬ 
culations.  This  method  is  still  widely  used  to  select 
boilers  and  chilleis;  however,  solar  energy  systems 
differ  from  conventional  systems  in  three  key  respects. 

First,  incident  solar  energy  is  not  continuously 
available  and  may  not  be  available  at  all  during  peak 
loads  at  night  or  in  cloudy  weather.  Thus,  since  solar 
(Hieigy  may  net  be  able  le  meet  ad  demands  ecwwtni 
cally,  provisions  must  be  made  for  auxiliary  energy 
supply  systems.  Sizing  of  solar  system  components  is 
therefore  based  on  annual  system  performance,  rather 
than  peak  load  calculations.  The  auxiliary  system  is 
still  sized  fur  peak  loads. 


^Interim  Performance  Criteria  for  Soiar  Heating  and 
Cooling  Systems  and  Dwellings  (National  Bureau  of  Standards, 
January  1, 1975). 
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Second,  holh  the  energy  demand  and  the  thermal 
energy  storage  mass  inlluence  the  temperature  of  tlie 
lliiid  entering  the  solar  collectors,  and  therefore 
inlluence  the  amount  of  incident  solar  energy  that  can 
he  collected  by  a  given  collector  array.  Consequently, 
seasonal  building  heating  and  ccnrling  load  variations 
strongly  influence  the  annual  performance  of  a  solar 
energy  system. 

Third,  solar  energy  is  typically  collected  and  deliv¬ 
ered  over  a  large  temperature  range.  This  inlluences 
the  performance  of  heating  coils  and  absorption  chil¬ 
lers.  Cor  these  reasmis,  conventional  “peak  load” 
design  methods  are  not  acceptable,  and  computer 
simulation  must  he  used. 

To  properly  .select  solar  system  components;  to 
establish  optimum  collector  area,  tank  volunK.  and 
collector  tilt  and  a/imuth  angle;  and  to  determine  the 
economic  feasibility  of  a  solar  heating  and  cooling 
system,  an  accurate  method  of  estimating  the  annual 
performance  of  a  given  solar  energy  system  must  be 
used.  .Such  a  method  is  provided  by  computer  simula¬ 
tion  techniques  which  use  hourly  weather  data  to 
predict  hourly  building  heating  and  cooling  loads, 
and  the  performance  of  candidate  solar  heating  and 
cooiit'g  systems.  The  computer  greatly  facilitates  the 
tedious  task  of  estimating  hourly  building  load  pro¬ 
files  and  performing  the  hourly  iterative  heat  balance 
calculations  required  to  estimate  annual  system  perfor¬ 
mance. 

It  is  recognized  that  load-predicting  and  solar  energy 
system  simulation  programs  are  not  yet  readily  avail¬ 
able.  Consequently,  parametric  computer  simulation 
studies  have  been  performed,  and  the  results  have  been 
used  to  develop  a  method  which  does  not  require  the 
further  aid  of  the  computer  to  estimate  solar  energy 
system  performance  for  buildings.  This  method  (de¬ 
scribed  in  Chapter  4)  was  developed  by  the  following 
procedure; 

(1)  Two  typical  Army  buildings  and  five  geogra¬ 
phical  sites  were  selected  for  study. 

(2)  Hourly  building  load  estimates  were  made 
for  each  building  at  each  site,  using  hourly  climato¬ 
logical  data  tapes  from  the  National  Climatic  Center, 
National  Oceanic  and  Atmospheric  Administration 
(NOAA),  Asheville,  NC,  and  the  National  Bureau  of 
Standards  Load  Determining  Program  (NBSLD). 
(Appendix  B  provides  a  general  desaiption  of  NBSLD.) 


Representative  weather  data  years  were  selected  for 
each  site. 

(.^)  Using  hourly  incident  solar  radiation  data  from 
the  NOAA  tapes  and  the  hourly  load  estimates  deter- 
ntined  above,  a  series  of  parametric  studies  was  per¬ 
formed  by  a  solar  sy.stem simulation  program  developed 
at  the  U  S.  Army  Construction  Hngineering  Research 
Uboratory  (CLRL).  These  parametric  studies  provided 
performance  estimates  of  the  various  solar  heating  and 
C(H>ling  systems. (.See  Appendix  C  for  a  detailed  descrip¬ 
tion  of  the  CERL  program.) 

Building  and  Syitam  Simulation 

The  study  investigated  two  structures  of  standard 
Army  design  which  are  typical  of  conventional  office 
and  dormitory  buildings:  a  modular  barracks  and  a  two- 
battalion  headquarters  and  classroom  building. 

Tire  barracks  (function  category  721 1 1)  is  a  three- 
story,  brick -faced  masonry  building  having  a  flat, 
built-up  roof  and  eight  .Tman  rooms  on  each  Hoor.  The 
72-man  modules  arc  typically  arranged  in  groups  of 
three  to  five  to  form  a  larger  barracks  complex.  (For 
this  study,  only  one  module  was  examined.)  Each 
module  contains  approximately  4(XX)  sq  ft  (372  m^) 
of  floor  space.  Ventilation  was  assumed  to  be  by  infil¬ 
tration  only. 

The  classroom  and  headquarters  facility  (function 
categoiy  61041)  is  a  single-story,  brick-faced  concrete 
block  structure  with  a  nearly  flat,  built-up  roof. 
Approximately  two-thirds  of  the  floor  space  contains 
windowless  classrooms.  The  remaining  one-third, 
which  IS  used  as  office  space,  has  20  exterior  windows. 
This  building  contains  approximately  12,000  sq  ft 
( 1 1 1  .S  m^ )  of  floor  space. 

Both  buildings  were  assumed  to  comply  with  the 
I  October  1972  DOD  Construction  Criteria  Manual 
4270.  I-M,  which  specified  a  unit  thermal  conductance 
minimum  of  .15  Btu/hr-sq  ft  (1.70  KJ/hr-m*)  for  walls, 
.10  Btu/hr-sq  ft  (1.1.3  KJ/hr-m^)  for  floors,  and  .05 
Btu/hr-sq  ft  (.57  KJ/hr-m*)  for  ceilings  and  roofs. 
The  buildings  were  assumed  to  be  heated  only  when 
the  space  temperature  was  68°F  (I9.8°C)  or  below, 
and  cooled  only  when  the  space  temperature  was  78°F 
(75  3°C)  or  above. 

The  buildings  described  above  were  analyzed  for 
five  sites  chosen  as  being  representative  of  the  various 
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climule  types  in  whieh  there  lire  liirtte  numbers  ol 
Army  buildings:  Hurt  Worth,  TX  (hot,  sunny  summers 
«nd  mild  winteTs),  l.ns  Ango1e>.  <  A  (high  im»rlirtimi 
jiid  relutively  unilorm  Icmpcratuies  thniughoul  the 
eeirl’  f'jiUunbus  Mfl  nsrld  wunl^ni  ku 

summers,  and  moderate  insolation):  Madison,  Wl 
(cold  winters,  fairly  short,  wami  summers);  and  Wash* 
ington,  IX'  (hot,  humid  summers,  too)  winters,  and 
moderate  insolation). 

I'lgtir''  I  shows  the  conliguration  of  the  modeled 
solar  heating  and  cooling  system.  Must  simulations 
were  performed  for  a  single-cover,  selective-surlace 
(.v^K-tos  Xiill'k'ieui  uuuilaiiwbi  wfu  nuule  ariilh  oclsei 
collectors  to  obtain  a  simple  conversi(..i  factor  to 
relate  the  performance  of  other  collector  types  to  the 


base  type.  (.See  the  “Universal  Curve  for  Estimating 
.Solar  Heating  and  CiMiling  Performance"  section, 
page  I S .)  The  thermal  energy  stmage  tank  was  assumed 
to  be  insulated  with  the  equivalent  of  approximately 

1  (n  ti  cni)  esC  poly:iHatluMK  UiMlIsiliioa  A  lUTidavl 
economy  cycle. was  assumed  in  the  calculations:  outside 
air  provided  all  cooling  when  its  temperature  was  less 
than  55"E  1 1 2.(>“C)  for  tlic  headquuirtrs,  ut  less  than 
65“F  ( IK.|“(')  for  the  barracks. 

Results  of  Simulation  Studies 

Table  I  sluiws  the  results  of  the  load  and  solar 
simulation  studies  for  the  barracks  module,  and  Table 

2  ihuws  re.uilu  tcur  tbe  he.«i4uiuxi;ni  budilatg  Valujer 
shown  for  optimum  collector  tilt  angle  were  obtained 
from  computer-derived  curves,  as  shown  in  f  igure  V 


Table  I 

Bamcki  Module 


Sin 

lORT  WORTH, 
TX 

(OLUMBIA, 

MO 

MADISON. 

Wl 

WASIHNOiON. 

D( 

LOS  ANOI  I.I  S, 
(A 

Year  uf  Weather  Data 

1955 

1965 

1961 

1954 

196.7 

I  alitude  of  Site  in  Degrees 

32 

39 

43 

39 

34 

Optimum  Collei'lur  Till  in  Degrees 

20 

26 

35 

25 

29 

Annual  Solar  Radiation  on  Col- 

6.54 

6.03 

6.23 

5.64 

6.68 

lector  al  Optimum  Tilt  (lU^  Htu/ 

(74.3) 

(68.5) 

(70.7) 

(64.0) 

(75.8) 

sq  ft)  ( 10’  ICJ/m*) 

Annual  llori/ontal  Solar  Radi- 

J.9 

5.2 

5.1 

4.9 

5.8 

alion  (to’  Btu/sq  fo  (lo’ 

(59.0) 

(58.0) 

(55.6) 

(65.0) 

Annual  Building  Heating  I.oad 

O.OOM 

0.40 

0.95 

0.16 

~0 

(10*  Blu/yr)(10'*  KJ/yr) 

(.008) 

(0.42) 

(1.00) 

(.17) 

(~0) 

Annual  Building  Cooling  I.oad 

4.16 

2.47 

1.36 

2.54 

t.76 

(10*  Btu/yr)(l0*  KJ/yr) 

(4.38) 

(  2.60) 

(  1.43) 

(-2.68) 

(  1.86) 

.Annual  Thermal  Knergy  Required 

6.4 

3.8 

2.1 

3.9 

2.7 

for  Absorption  Cooling  With 

(6.7) 

(4.0) 

(2.2) 

(4.1) 

(2.8) 

COP=0.6S  (10®  Btu/yr)(10" 

KJ/yr) 

Total  Annual  Knergy  Required 

6.4 

4.2 

3.03 

4.15 

2.70 

(10*  Btu/yr)(10*  KJ/yr) 

(6.7) 

(4.4) 

(3.1) 

(4.3) 

(2.8) 

Pealc  Heating  Load  ( lo’  Btu/hr) 

0.32 

0.89 

1.3 

0.065 

~0 

(10*  KJ/hr) 

(.34) 

(.94) 

(1.4) 

(.068) 

(0) 

Peak  Cooling  Load  ( 1 0*  Blu/hr) 

1.42 

-1.30 

-1.20 

1.36 

1.28 

(10*  KJ/hr) 

(1.44) 

(-1.36) 

(-1.27) 

(-1.43) 

(-1.35) 
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Figure  3.  Effect  of  collector  tilt  angle  on  fraction  of  load  met  by  solar  energy. 


Tabic  2 

Haadqurtcrt  BuiMiiig 


SITE 

l  ORT  WORTH, 
TX 

COLUMBIA. 

MO 

MADISON. 

Wl 

WASHINGTON, 

IX 

LOS  ANGELES, 
CA 

Optimum  Colleclor  Till  (De(reeil 

22 

28 

28 

27 

31 

Annuiil  Subr  Radialion  on  C'ollectoi 

6.55 

6.05 

6.18 

5.67 

6.69 

at  (^timum  Till  (10^  Blu/sq  ft) 

(74.4) 

(68.7) 

(70,2) 

(64.4) 

(76.0) 

Annual  Building  Heating  Load 

0.4 

1.0 

1.8 

0.8 

0.4 

(10*  Blu/yr)(IO*  KJ/yr) 

(4) 

(1.1) 

(1.9) 

(.8) 

(.4) 

Annual  Building  Cooling  Load 

-7.7 

5.3 

3.5 

5.0 

-5.7 

(IO*Blu/yr)(IO*  KJ/yr) 

(  8.1) 

(5.6) 

(-3.7) 

(  5.3) 

(-6.0) 

Annual  Thermal  Energy  Required 

11.9 

8.) 

5.4 

7.7 

8.8 

fur  Absorption  Cooling  With 

(12.5) 

(8.5) 

(5.7) 

(8.1) 

(9.3) 

OOP -0.65  (10*  Btu/yr) 

(10*  KJ/yr) 

Total  Annual  Energy  Required 

12.34 

9.08 

7.21 

8.52 

9.19 

(10*  Btu/yr)  (10*  KJ/yr) 

(13.0) 

(9.5) 

(7.6) 

(9.0) 

(9.7) 

Peak  Heating  Load  ( 10*  Btu/hr) 

— 

0.58 

0.89 

0.24 

_ 

(10*  KJ/hr) 

(.61) 

(.94) 

(.25) 

Peak  Cooling  Load  (10*  Btu/hr) 

-3.0 

-2.9 

2.8 

2.9 

-2.8 

(10*KJ/hi)  -- 

(3.2) 

(-3.0) 

(2.9) 

(-3.0) 

(-2.9) 

I  hose  curves  show  iliut  the  tilt  uiigle  is  not  critical;  a 
deviation  ol  15"  Iroin  optimum  has  little  cITecl  on 
IH'rlormancc.  It  is  iiitcresliii^  to  note  that  when  the 
collector  has  an  appropriate  tilt,  latitude  Is  not  a 
domitialuii!  lactor  ariectint!  annual  incident  radiation. 

At  all  sites,  even  the  coldest,  the  cooling  htad  ex¬ 
ceeds  the  heatin(>  load.  The  line  in  Tables  I  and  2  labeled 
“Annual  Thermal  Knerny  Required  lor  Absorption 
C’cHiling  With  COP  =  O.h.S”  is  the  thermal  energy  required 
by  the  absorption  cooler  to  meet  the  building  cooling 
load.  Adding  the  healing  load  to  this  value  gives  the 
total  annual  thermal  energy  required. 

Figures  4-N  are  plots  of  the  solar  perlormance  curves 
lor  the  barracks  module  of  each  sample  site.  The  ordi¬ 
nate  on  each  graph  is  the  tract  ion  ot  the  total  annual 
thermal  load  supplied  by  solar  ertergy  tor  the  particular 
site.  The  abscissa  is  the  total  ellective  collector  area. 
Kach  graph  has  lour  curves  corresponding  to  dilferent 
ratios  of  storage  water  mass  to  collector  area 
where  is  in  pounds  of  water  and  A^.  is  in  square 
feel.  (One  lb  of  water  per  s<|uare  fool  of  collector  area 
corresponds  to  4.‘)  kg  of  water  per  square  meter  of 
collector.)  These  curves  are  representative  of  the  type 
of  performance  curve  that  must  be  constructed  to  per¬ 
mit  determination  of  the  economic  feasibility  of  solar 
heating  and  cooling. 

Figure  9  compares  the  solar-supplied  energy  for 
different  sites  for  an  arbitrarily  selected  storage  capac¬ 
ity  of  .41  lb  water  per  square  f(K>t  of  collector  area 
(IS2  kg/m^|.  It  can  be  seen  that  for  a  system  of  given 
si/e,  the  anxrunt  of  solar-supplied  energy  ranges 
greatly  for  the  different  sites,  with  Fort  Worth  having 
the  greatest  amount  and  los  Angeles  the  smaOest. 
These  differences  are  largely  due  to  the  differences  in 
total  annual  thermal  energy  demanded  by  the  building 
at  each  site  and  indicate  the  critical  impact  of  building 
energy  use  on  the  overall  performance  of  the  solar 
energy  system.  Figure  10  is  an  example  of  how  the 
fraction  of  load  provided  by  solar  energy  depends  on 
the  thermal  storage  mass,  with  collector  area  as  a  para¬ 
meter.  These  curves  show  a  big  improvement  for 
increasing  storage  from  K  to  16  lb  water/sq  ft  (38  to 
.304  Kg/m^)  and  only  minor  improvement  for  increasing 
storage  from  31  to  62  Ib/sq  ft  (152  to  304  Kg/m^). 

Figures  II  to  IS  are  performance  curves  fur  the 
headquarters  building  for  the  five  sample  sites.  The 
format  is  the  same  as  the  curves  for  the  barracks 
module.  Figure  16  is  a  comparison  of  the  sites. 


Universal  Curve  for  Estimatinfl  Solar 
Heating  and  Cooling  Parformanea 

Figure  17  is  a  dimensionless  comparison  of  both 
buildings  at  all  sites.  The  abscissa  is  the  incident  annual 
radiation  on  the  collector  with  optimum  tilt  (the 
product  of  collector  area  and  radiation  density)  divided 
by  the  total  annual  energy  requirements.  The  ordinate 
is  the  fiaction  of  this  annual  energy  requirement  met 
by  solai  energy.  This  plot  is  for  a  storage  of  1 6  lb  water 
per  square  foot  of  collector  area.  The  simulation  data 
points  fall  remarkably  close  to  the  single  representative 
curve  shown.  Thus,  Figure  17  can  provide  a  reasonably 
good  approximation  of  solar  heating  and  cooling  sys¬ 
tem  performance  fur  all  buildings  at  all  sites.  Figures 
18  and  19  are  similar  plots  with  representative  curves 
for  larger  storage  sizes.  Figure  20  shows  the  three 
representative  curves.  The  constants  were  obtained  by 
a  computer  determination  of  the  best  “least  square” 
fit  of  tlic  simulated  data. 

Each  curve  is  described  by  the  general  equation 

p  *  B(r- 0.082  r*)  (Eql) 

where  p  is  the  fraction  of  the  annual  load  supplied  by 
solar  energy,  and  r  is  a  dimensionless  parameter  defined 
by 


where; 

is  the  annual  radiation  energy  per  unit  area  on 
the  tilted  collector, 

A  is  the  collector  area, 

c 

M  is  the  multiplying  factor  of  Table  3 ,  and 
is  the  total  annual  energy  requirements. 

The  constant  B,  whose  values  are  shown  in  Figure 
20.  depends  on  the  chosen  storage  size. 

If  values  for  H.^,  Q^,  B,  and  M  are  known,  Eq  I 
reduces  to  p  as  a  function  of  the  collector  area,  A^, 
alone-the  ume  form  as  the  performance  curves 
shown  earlier.  This  equation  can  then  be  used  for 
economic  feasibility  calculations. 

The  annual  radiation  density  on  the  tilted  collector 
surface  may  be  obtained  by  using  the  empirically 


Figure  4.  Performance  curve  for  barracks  at  Fort  Worth,  TX,  with  total  annual  energy  requirement  of 
6.4  X  10*  Btu  (6.7  X  10"  joules). 
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Figure  5.  Performance  curves  for  barracks  at  Washington,  DC,  with  total  annual  energy  requirements  of 
4.05  X  10*  Btu  (4.27  X  lO' '  joules). 
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Figure  6.  Performance  curves  for  barracks  at  Columbia,  MO.  with  total  annual  energy  requirements  of  4.2  X  10* 
Btu(4.4  X  10* '  joules). 
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Figure  7.  Performance  curves  for  barracks  at  Madiion,  Wi,  with  total  annual  energy  requirements  of  3.03  X  I  O' 
Btu(3.20X  10"  joules). 
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Figure  8.  Performance  curve  for  barracks  at  Los  Angeles.  CA,  with  total  annual  energy  requirements  of  2.7  X  10 
Btu(2.8SX  10"  joulei). 
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Figure  9.  Compiriton  of  performance  curves  for  barracks  at  all  sites. 
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Figure  10.  hrfecl  ol  thermal  storage  capacity  on  performance. 
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Figure  II.  Performance  curves  for  headquarters  building  at  Fort  Worth, TX,  with  total  annual  energy  requirements 
on2.34X  10*  Btu  (13.02  X  lO"  joules). 
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Figure  12.  Performance  curves  for  headquarters  building  at  Columbia,  MO,  with  total  annual  energy  requirements 
ofd.OSX  IO''Btu(9.58X  lO" joules). 
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F%ui«  13.  Performance  curves  for  headquarters  building  at  Madison,  WI,  with  total  energy  requirements  of  7.21  X 
10*  Btu  (7.61  X  10' '  joules). 
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Figure  14.  Pertormance  curves  for  headquarters  building  at  Washington,  DC,  with  total  annual  energy  requirements 
of  8. 52  X  10"  Btu(<i.0X  I O"  joules). 
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Figure  IS.  Performance  curves  for  headquarters  building  at  Los  Angeles,  CA,  with  total  annual  energy  requirements 
of  9.19  X  10*  Btu(9.7X  10* ‘joules). 
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Figure  16.  Comparison  of  performance  curves  for  headquarters  building  at  all  sites  with  storage  =  31  lb  water.'sq  ft 
(152kg/m^). 
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Figure  17.  Universal  performance  curve  for  both  buildings  at  all  sites.  Storage  *  16  lb  water/sq  ft  ( 76  kg/m^ ). 
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Figure  18.  Universal  performance  curve  for  both  buildings  at  all  sites.  Storage  31  lb  water/sq  ft  (152  kg/m^). 
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ANNUAL  INODENT  RADIATION/ANNUAL  ENERGY  REQUIREMENTS,  r 
Figure  19.  Universal  performance  curve  for  both  buildings  at  all  sites.  Storage  =  62  lb  water/sq  ft  (304  kg/m^ ). 
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Figure  20.  Approximate  universal  performance  curves  for  three  storage  volume-to-collector  area  ratios  obtained 
from  the  analytic  expression  shown. 
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derived  Ibrmula 


M 


ciis(0|  X) 


(N .») 


where 

M  is  Ihe  atiiuiul  rudiation  density  on  a  liori/onlal 
surface  and  is  available  for  most  regions  in  llic 
United  States,  and 

0^  is  the  latitude  of  the  chosen  site. 

A  map  of  mean  daily  horizontal  radiation,  taken 
from  the  Climatic  Atlas  of  the  United  S'ates  (U.S. 
Dept  of  Commerce)  is  shown  in  Figure  21.  The  radia¬ 
tion  density  is  given  in  Langleys/day.  The  annual 
radiation  density  in  Btu/sq  ft  is  obtained  by  multi¬ 
plying  the  value  from  the  map  by  I  ..^4  X  10^ ;  to  obtain 
11^1  in  joules/m^ ,  the  factor  is  I  ..S3  X  10^. 

The  problem  of  determining  the  annual  heating  and 
tooling  loadtO^*  addressed  on  p  27. 

Since  all  the  simulations  were  made  with  a  collector 
of  particular  design  (single  cover  with  selective  coating; 
a*0.90;  e=O.IO,  n=l),  a  comparison  is  made  in  Table 
3  for  collectors  with  other  properties.  In  Table  .3  a  is 
the  absorptivity  of  visible  radiation  of  the  absorbing 
plate,  and  e  is  the  infrared  emissivity  of  the  plate. 
Values  of  a  and  e  are  available  from  the  collector 
manufacturer.  The  values  in  the  table  are  the  multiply¬ 
ing  factor,  M,  used  in  bq  2  to  obtain  performance 
curves. 

By  using  the  curves  of  Figure  20  in  conjunction 
with  the  adjustment  factors  for  collector  type,  build¬ 
ings  and  site-speciHc  curves  can  be  constructed  which 
show  the  fraction  of  annual  load  met  from  solar  energy 
as  a  function  of  collector  area.  Separate  curves  of  this 
type  can  be  drawn  for  different  storage  mass  to  collec¬ 
tor  area  ratios,  as  illustrated  by  Figures  1  through  8 
and  II  through  IS,  and  for  different  collector  types. 
Such  curves  permit  determination  of  economically 
optimum  designs,  as  will  be  described  in  Chapter  4. 

It  should  be  noted  that  for  buildings  with  dominant 
cooling  loads,  if  the  annual  heating  and  cooling  load 
of  the  building  and  the  annual  radiation  density  of  the 
site  are  known,  the  previously  described  performance 
estimating  method  can  be  applied  regardless  of  the 
type  of  building  or  site.  However,  the  method  applies 
only  to  systems  designed  to  meet  both  heating  and 


cooling  loads.  The  curves  do  not  apply  to  heating  only 
applications  or  to  applications  which  have  small  cool¬ 
ing  loads.  .Such  curves  are  being  prepared  and  will  be 
available  in  the  near  future. 

An  example  will  illustrate  how  the  universal  curves 
t)f  Figure  20  may  be  used  to  determine  the  perfor¬ 
mance  of  a  solar  installation.  Consider  a  building  at 
Nashville, TN,  which  is  known  to  have  an  annual  heating 
load  of  I  X  lO"  Btu  and  an  annual  cooling  load  of  4 
X  10*  Btu.  The  thermal  requirement  for  the  cooling  is 
then  4  X  I0*/0.6S  =  6.1  S  X  10*,  where  0.6.S  is  the 
coefficient  of  performance  of  the  absorption  cooler. 
The  total  annual  energy  requirement,  Q^,  is  then  I  X 
10*  +6. US  X  10*  =7.15  X  10*  Btu. 

From  the  map  for  solar  radiation  (Figure  21),  we 
find  that  for  Nashville,  the  horizontal  radiation  density 
is  3SS  Langleys/day,  and  the  latitude  is  about  32°. 
is  then  equal  to  355  X  1.34  X  10^  =  4.8  X  10*  Btu/sq  ft 


Table  3 

Collector  Area  Multiplying  Factor,  M 


COLLECTOR  AREA  MULTIPLYING  FACTORS 
FOR  DIFFERENT  COLLECTOR  DESIGNS 


0.96 

0.94 

0.90 

e 

0.96 

0.30 

0.10 

N 

1 

1.55 

1.09 

1 

2 

1.09 

0.97 

0.93 

ABSORPTIVfTY 
£•  EMISSIVITY 

N  •  NUMBER  OF  GLASS  COVERS 
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DAILY  SOLAR  RADIATION  (Langleys).  ANNUAL 


F|psc  21.  Solar  radiation  map  of  U.S.  From  Qimatic  Atlas  of  the  United  States  (U.S.  Department  of  Commerce, 
1%8). 


annually.  Using  Ivq  .V  llie  rudialion  ilensily  on  llic 
lilted  collector,  is  found. 

4.8  X  10* 

H  >  -  =  5.2.‘i  X  10*  Btu/su  ft 

®  cos  (.12  8) 

From  Fq  2,  assuming  M=l  (single  cover,  selective 
coating),  we  determine 

.S.2.SX  I0*A^, 

r  »  - ^  =  7..14X  10 

7.I5X  10* 

Using  this  relation,  Fq  I,  and  Figure  20,  and  selecting 
a  storage  of  31  Ih/sq  ft  ( I.SI  kg/m^)  the  frr-lion  of  the 
energy  requirements  supplied  hy  solar  gives  energy: 

0  =  0.309  I7.34X  10  “A^  0.08  2  (7..34X  IO  -’A^)^] 

=  2.27X  I0-*  (A  6.02XI0  *A^) 

Table  4  shows  values  of  p  for  a  range  of  collector  sizes. 


Table  4 

Energy  Requirements  Mel  by  Solar  Energy 


(«|  ft)  500 

1000 

1500 

2000 

2500 

3000 

3500 

fi  0.1 1 

0.21 

0.31 

0.40 

0.48 

0.56 

0.63 

The  table  contains  the  performance  data  needed  to 
make  an  economic  feasibility  assessment. 


Determining  Anmiel  Heating  and 
Cooling  Energy  Load 

It  is  not  the  objective  of  this  report  to  provide  a 
method  for  precise  determination  of  annual  building 
heating  and  cooling  energy  requirements;  however, 
the  problem  must  be  briefly  addressed  because  of  its 
impact  on  solar  energy  system  economic  feasibility 
assessment. 

For  sites  and  buildings  other  than  those  for  which 
annual  energy  requirements  are  given  in  the  “Results 
of  Simulation  Studies"  section,  p  13,  the  annual 
building  heating  and  cooling  energy  requirements  can 
best  be  estimated  from  measured  consumption  data  for 


identical  or  similar  buildings  near  the  building  being 
considered.  When  this  is  not  possible,  hourly  load- 
predicting  conipulcr  programs  which  use  actual  weath¬ 
er  data  are  the  next  best  choice.  The  NBSl.D  Program 
(Appendix  B)  and  NASA’s  "NFCAP"  Program  provide 
good  results.  If  other  programs  requiring  the  input  of 
peak  heating  and  c(K)ling  loads  arc  used,  these  loads 
should  be  calculated  with  response  factor  methods 
outlined  in  the  ASHRAB  Handbook  of  Fundamentals, 
1973  edition. 

When  the  techniques  described  above  arc  not  avail¬ 
able,  a  rough  approximation  of  the  heating  and  cooling 
load  can  be  obtained  by  scaling  the  results  presented 
in  the  “‘Results  of  Simulation  Studies”  section.  The 
first  step  in  this  scaling  procedure  is  to  adjust  for  sites 
not  l(x;atcd  at  or  near  the  five  studied.  Various  ap¬ 
proaches  have  been  suggested  for  adjusting  the  annual 
heating  and  cooling  load  to  account  for  climate  differ¬ 
ences.  Most  correlations,  however,  have  been  poor. 
In  particular,  peak  healing  and  crxrling  loads  have  little 
or  no  relationship  to  geographical  differences  in  annual 
energy  consumption.  Degree  days,  particularly  “cool¬ 
ing  degree  days"  are  also  poor  adjustment  factors. 
Therefore,  intuitive  adjustments  based  on  information 
discussed  in  the  “Results  of  Simulation  Studies" 
section  appear  to  be  the  best  recourse.  Interpolation, 
weighted  by  judgment,  between  annual  loads  for  the 
two  sites  closest  in  climate  to  the  location  being  con¬ 
sidered  appears  to  be  the  best  approach.  For  example, 
for  estimating  the  annual  load  of  a  barracks  in  central 
Illinois,  interpolation  between  the  loads  of  Columbia, 
MO  and  Madison,  Wl,  should  produce  a  reasonable  re¬ 
sult. 

When  site  adjustments  have  been  made,  adjustments 
for  differences  in  building  size  or  configuration  must 
be  accomplished.  Note  that  the  buildings  selected  for 
the  computer  simulation  study  were  chosen  because 
they  typify  two  of  the  major  functional  building  types 
found  at  Army  posts,  (i.e.,  troop  housing  and  admini¬ 
strative  buildings)  and  are  representative,  in  terms  of 
their  thermal  responses,  of  most  concrete  masonry 
and  brick  construction.  Consequently,  the  results  of 
the  energy  requirement  calculations  for  two  buildings 
at  five  sites  presented  in  the  “Results  of  Simulation 
Studies"  section,  (p  13),  may  be  legitimately  extended 
to  other  masonry  and  brick  construction  by  applying 
a  rough  scaling  factor  based  on  the  ratio  of  the  sum  of 
roof  and  exterior  wall  area  of  a  candidate  building  to 
that  of  one  of  the  buildings  studied.  This  scaling  should 
not  be  applied,  however,  to  lightweight  or  uninsulated 
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builJiiigs  ot  lo  hiiildintts  with  1ur^c  ({lass  areas,  siiu'o 
the  ihoriual  properties  ol  these  buildings  diHer  largely 
Initri  those  stitdied. 

Note  that  regardless  of  the  iiK'thod  used  to  deter- 
iiiMie  lire  annual  building  energy  rei|uiretnents,  the 
value  required  to  implement  the  procedures  discussed 
rti  the  ‘‘Dnivcrsal  Curve’’  section  is  not  the  sum  of  the 
absolute  values  of  the  annual  hitilding  heating  load  and 
the  annual  cooling  load;  it  is  rather  the  amount  of 
etierg>  reijuired  by  the  healing  and  cooling  systems 
to  meet  the  building  load  (note  the  difference  between 
cooling  loads  and  energy  required  l«)r  cooling  in 
Tables  I  and  2).  This  distinction  is  critically  important, 
not  only  because  of  the  cocfncicnl  of  t<crfurmance 
ix'nally  associated  with  absorption  ccxtling,  but  also 
because  large  energy  consumption  differences  can 
occur  between  different  types  of  fan  systems  when 
meeting  identical  building  loads  (compare,  fur  exam¬ 
ple,  two-pipe  fan  coil  versus  terminal  reheat  systems). 

This  leads  to  another  general  rule  which  must  be 
considered  in  the  application  of  solar  energy.  It  is 
currently  more  cost-effective,  in  almost  every  case, 
to  apply  strict  energy  conservation  techniques  to  a 
building  before  considering  solar  energy  as  an  alter¬ 
native  means  of  heating  and  cooling  it.  The  applica¬ 
tion  of  solar  energy  should  therefore  be  considered  in 
its  proper  perspective  as  one  of  many  mutually  com¬ 
patible  energy  conservation  options. 


4  ECONOMIC  FEASIBILITY 

ASSESSMENT  AND  SELECTION 
OF  OPTIMUM  SYSTEM 

The  method  described  in  this  chapter  for  defining 
the  economically  optimum  solar  energy  system  relies 
both  on  the  building’s  annual  energy  consumption 
(p  27)  and  on  the  solar  energy  system  performance 
estimates  (pl5)  described  in  Chapter  3.  The  proce¬ 
dure  is  an  adaptation  of  that  suggested  by  Butz.et  at.* 
and  leads  to  the  selection  of  a  collector  area  and  lank 
volume  which  will  maximize  the  “trade-off'  between 
solar  energy  system  costs  and  fuel  cost  savings.  The 
method  relies  on  a  simple  graphical  comparison  of  the 
life-cycle  (present  worth)  cost  difference  between  a 


*L.  W.  Butz,  et  at..  Use  of  Solar  Unergy  for  ResUential 
Heating  and  Cooling,  M.  S.  Thesis  in  Mechanical  Enfincering 
(Univeriityof  Wiicontin,  1973). 


solar  energy  system  and  a  conventional  system.  Appen¬ 
dix  A  provides  an  example  application  of  the  method 
to  a  hattalion  headi|uarlors  and  classroom  building  at 
l-ort  Hood,  rX. 

The  first  task  is  estimation  of  the  lile-eyele  cost  ol 
the  solar  energy  system  components.  Because  this 
method  is  comparative,  only  the  costs  of  components 
that  are  not  normally  |iart  of  a  conventional  healing 
and  cooling  system  should  be  considered;  thus,  the 
cost  of  the  building's  air-handling  system  would  not 
He  considered,  but  the  difierence  between  the  installa¬ 
tion  cost  of  a  mine  expensive  absorption  chiller  and  of 
a  less  expensive  centrifugal  chiller  should  be  charged 
to  the  solar  energy  system.  Note  that  for  a  solar 
heating  and  cooling  system,  certain  cost  elements  vary 
according  to  the  size  of  solar  energy  systems,  and 
others  are  relatively  fixed,  regardless  of  the  collector 
area  or  tank  volume  chosen. Collector  and  storage  tank 
costs  arc  obvious  examples  of  systcm-sizi. -dependent 
items;  other  examples  include  heat  exchanger  costs, 
and  certain  pump  and  piping  costs.  Hie  additional 
control  system  cost  associated  with  a  solar  energy 
system  is  an  example  of  a  cost  difference  that  is  largely 
independent  of  collector  area.  The  cost  difference 
associated  with  the  purchase  and  installation  of  an 
absorption  chiller  is  also  relatively  independent  of  solar 
collector  area,  since  for  all  but  the  smallest  solar 
collector  areas,  selection  of  an  appropriate  absorber  is 
dictated  by  the  peak  building  cooling  load. 

When  these  fixed  and  variable  cost  items  are  separ¬ 
ated,  it  is  a  relatively  straightforward  matter  to  esti¬ 
mate  the  life-cycle  cost  difference  between  solar  and 
conventional  heating  and  cooling  systems  for  various 
collector  area/storage  capacity  combinations.  Figure 
22  illustrates  how  the  capital  component  of  the  life- 
cycle  cost  increases  with  collector  area. 

The  second  step  of  the  comparison  procedure  is 
establishment  of  the  life-cycle  fuel  costs  for  conven¬ 
tional  and  solar  energy  systems  of  varying  sizes. 

For  the  conventional  system,  annual  fuel  costs  can 
be  estimated  as  follows; 

(1)  Determine  the  annual  heating  and  cooling 
energy  loads  for  conventional  equipment,  using  the 
procedures  discussed  in  Chapter  3,  p  27. 

(2)  Convert  the  cooling  energy  load  to  kilowatt 
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COLLECTOR  AREA  (FOR  FIXED  STORAGE  TO  COLLECTOR 
AREA  RATIO) 


Figure  22.  Life-cycic  cost  difference  vs.  collector  area. 


hours,  and  multiply  by  the  local  electrical  rate  to  ob¬ 
tain  the  cost  of  cooling  energy. 

(3)  Divide  the  input  heat  energy  required  for 
heating  by  the  heating  value  of  the  fuel  used,  and 
multiply  by  the  unit  fuel  price  to  determine  the  cost 
of  heating. 

(4)  Add  the  results  of  (2)  and  (3)  to  obtain  the 
total  annual  fuel  cost. 

(5)  Convert  the  annual  fuel  cost  to  life-cycle 
(present  worth)  cost. 

The  amount  of  auxiliary  fuel  or  electrical  energy 
required  annually  for  heating  and  cooling  can  be 
determined  for  various  solar  energy  system  sizes  by 
referring  to  the  site  and  building-specific  performance 
curves  generated  by  methods  outlined  in  Chapter  3. 
p  IS.  To  simplify  the  economic  analysis,  it  is  generally 
assumed  that  all  thermal  energy  requirements  are  met 
by  a  single  fuel  source  (gas,  oil,  or  electricity).  The 
procedure  for  determining  life-cycle  costs  for  various 


solar  energy  systems  is: 

(1)  Determine  the  total  auxiliary  annual  energy 

requirements  for  a  given  system,  Q^.,  by  multiplying 
the  total  energy  requirements,  Q^,  by  the  fraction  of 
the  total  requirements  not  met  by  solar  energy,  ie., 
0,  =  0l(1  P). 

(2)  Determine  annual  energy  cost  by  dividing  the 
annual  energy  requirements  by  the  heating  value  of  the 
fuel,  and  multiplying  this  figure  by  the  appropriate 
unit  fuel  price. 

(3)  Convert  annual  fuel  cost  to  life-cycle  (present 
worth)  cost. 

Figure  23  is  a  hypothetical  plot  of  life-cycle  fuel 
costs  versus  collector  area  generated  by  repeating 
steps  1  through  3  for  different  collector  areas. 

The  total  life-cycle  cost  difference  between  conven¬ 
tional  and  solar  heating  and  cooling  systems  can  now 


29 


I 


—  ■""■SHiii 


COLLECTOR  AREA  (FOR  FIXED  STORAGE  TO  COLLECTOR 

AREA  RATIO) 


Figure  23.  Life<ycle  fuel  costs  vs.  collector  area. 


be  determined  for  various  collector  areas.  For  a  given 
collector  area,  this  difference  is  simply  the  life<ycle 
capital  cost  difference  of  the  solar  energy  system 
(Figure  22)  plus  the  life-cycle  fuel  cost  for  the  same 
system  (Figure  23)  minus  the  life-cycle  fuel  cost  of  a 
conventional  system.  Figure  24  is  a  hypothetical  plot 
of  this  cost  difference  for  varying  collector  areas. 

When  the  cost  difference  for  the  solar  energy  sys¬ 
tem  being  considered  has  been  plotted  (as  in  the  form 
of  Figure  24),  the  optimum  collector  area  becomes 
obvious.  To  determine  the  optimum  solar  energy  sys¬ 
tem,  separate  building  and  site-speciflc  curves  should 
be  drawn  for  various  collector  types  and  storage- 
mass-to-coUector  area  ratios.  Note  that  Figure  24 
concisely  summarizes  the  economic  feasibility  of  solar 
heating  and  cooling  systems.  The  points  on  the  curve 
having  a  positive  cost  difference  indicate  that  solar 
energy  is  not  cost-effective.  Only  when  the  curve  dips 
below  the  origin  can  a  solar  energy  system  be  economi¬ 
cally  justified.  Based  on  present  solar  energy  system 
costs,  many  applications  may  not  have  such  a  negative 
cost  difference. 


5  CONCLUSIONS  AND 
RECOMMENDATIONS 

Conduiions 

Solar  energy  can  meet  much  of  the  heating  and 
cooling  loads  of  Army  buildings  at  most  locations 
within  the  United  States,  and  thus  can  potentially 
lower  a  facility’s  demand  for  conventional  fuels  signif¬ 
icantly.  Technically  feasible  solar  energy  systems  can 
be  installed  at  Army  facilities  with  currently  available 
commercial  components. 

The  economic  viability  of  solar  energy  systems 
varies  greatly  with  location  and  building  structure; 
thus,  each  application  must  be  analyzed  individually. 
Study  results  indicate  that  the  single  dimensionless 
curves  shown  in  Chapter  3  are  reasonable  approxima¬ 
tions  of  volar  heating  and  cooling  system  performances 
for  all  buildings  and  sites  studied.  Therefore,  for 
feasibility  analysis,  solar  system  performance  can  be 
estimated  without  using  computer-aided  techniques. 
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The  meihuds  described  herein  will  provide  a  first 
estimate  ut' economic  feasihiliiy;  if  this  first  estimate 
is  promising,  a  mure  specific  and  detailed  study  can 
provide  more  precise  design  information. 

An  unavoidable  conclusion,  however,  is  lliat  the 
design  of  solar  energy  heating  and  c(H)ling  systems  is 
loo  complex  for  handbook  methods;  computer  simula¬ 
tions  of  the  form  presented  in  this  report  are  therefore 
required  to  obtain  adequate  design  solutions. 

Racommandationi 

It  is  recommended  that  solar  energy  be  considered 
not  as  a  panacea,  but  as  one  of  many  viable  energy 
conservation  alternatives  which  must  be  applied 
collectively  to  reduce  the  Army’s  and  the  nation’s 
wasteful  use  of  dwindling  energy  resources. 


cation,  Mt’lhnds  a)  Testing  /hr  Rating  Solar  Collectors 
Based  on  Thermal  Performance^  should  be  applied 
whenever  the  manufacturer's  performance  claims  are 
in  question. 

.Solar  healing  and  cooling  systems  should  generally 
be  subject  to  the  same  reliability  standards  as  conven¬ 
tional  healing  and  cooling  systems.  The  NBS  publica¬ 
tion,  Interim  Performance  Criteria  for  Solar  Heating 
and  Cooling  Systems  and  Dwellings,'’  provides  add¬ 
itional  information. 

The  bnergy  Research  and  Development  Administra¬ 
tion,  in  cooperation  with  the  Department  of  Housing 
and  Urban  Development,  has  prepared  a  directory  of 
manufacturers  of  solar  heating  and  cooling  system 
components,^  which  may  be  used  as  a  source  for 
current  solar  energy  system  component  cost  data. 


It  is  recommended  that  the  method  described  herein 
be  field-tested  at  a  Corps  District  Office  to  validate 
its  ease  of  use  and  applicability  to  designers.  It  is  also 
recommended  that  the  information  and  methods 
described  in  this  report  be  applied  to  all  new  construc¬ 
tion  and  to  most  existing  well-insulated  buildings  to 
determine  the  economic  feasibility  of  using  solar 
energy  for  heating  and  cooling. 

When  comparing  various  types  or  brands  of  solar 
collectors,  the  procedures  described  in  the  NBS  publi- 


^J.  S.  Hill  and  T.  Kusuda,  Methods  of  Testing  for  Rating 
Solar  Collei  tors  Rased  on  Thermal  Rerformamv,  MtSIR-74- 
635  (Thermal  Engineering  Systems  Section,  Center  for  Build¬ 
ing  Technology,  December  1974). 

'’Interim  Perform  jnce  Criteria  for  Solar  Heating  and  Cool¬ 
ing  Systems  and  'swellings  (National  Bureau  of  Standards, 
January  1, 1975). 

^Catalog  on  Solar  Knergy  Heating  and  Cooling  Products 
(l'.RI)A.  1975).  (Available  from  Technical  Information  Center. 
P.O  Box  62,  Oak  Ridge,  TN  37830.) 


Figure  24.  Net  life-cycle  cost  difference  vs.  collector  area. 


APPENDIX  A; 

EXAMPLE  OF  ECONOMIC 
FEASIBILITY  CALCULATIONS 

I  he  general  approach  to  determining  economic 
Icjsihility  was  described  in  (  haptcr  4,  p  2b,  As  an 
illustration,  a  spccil'ic  building  the  headquarters  and 
classroom  at  a  specific  site,  hort  Hood,  TX  has  been 
chosen  arbitrarily.  In  order  to  make  an  economic 
assessnKMit,  an  appropriate  performance  curve, p(A^,), 
must  be  selected.  Assuming  that  a  storage  capacity  of 
.fl  lb  for  each  sq  ft  (IS  I  kg/m^)of  collector  is  reason¬ 
able  lor  the  first  estimate,  we  find  from  Figure  20  that 

p  =  O..W)(r  0.082  r^). 

The  paranKti.M, depends  on  the  annual  energy  require¬ 
ments,  Q|  ,  ol  the  building  and  the  radiation  density. 
H^.  Since  it  is  to  be  expected  that  climatic  conditions 
at  Fort  Hood  arc  quite  similar  to  those  at  Fort  Wurth, 
the  Fort  Worth  value  from  Table  2,  Oj^  =  1.23  X  lO’ 
Btu,  can  be  used.  From  the  solar  radiation  map  (Figure 
21 1,  a  daily  horizontal  radiation  of  440  Langleys,  and 
a  latitude,  0^,  of  32°  is  estimated.  To  obtain  the 
annual  horizontal  radiation  density,  the  daily 
value  in  Langleys  is  multiplied  by  1.34  X  lO’  Btu/sq 
ft .  The  value  fur  H^  is  then 

H^ 

H.  - - - - -  6.46  X  10*  Btu/sq  ft. 

®  cos(0j^-8) 

This  leads  to  r  =  HgA^VO^  =  5.25  X  10  ^  X  A^.  Substi- 
tuting  this  into  the  equation  for  p  gives  the  perfor* 
inance  curve, 

p=  1.62  X  10'*(A^-4.3X  10  *A^^)  (Eq  Al) 

which  will  be  used  in  the  following  cost  analysis. 

Before  examining  the  costs  involved,  some  discus¬ 
sion  of  the  cost  data  included  in  this  example  is  re¬ 
quired.  Cost  figures  for  capital  equipment  were  gener¬ 
ally  obtained  from  equipment  manufacturers.  Other 
cost  data,  such  as  controls,  piping,  and  labor,  were 
based  on  the  authors’  experience  and  the  use  of  stan¬ 
dard  estimating  guides.  Obviously,  these  cost  figures 
will  vary  with  location,  manufacturer,  and  time;  thus, 
although  these  figures  are  presently  representative, 
they  should  not  be  taken  as  guidelines  for  economic 
studies.  In  addition,  since  this  wu  an  example  to 


demonstrate  the  method  ol  performing  the  economic 
analysis,  no  attempt  was  made  to  obtain  hard  economic 
data;  standard  textbook  present-worth  costing  methods 
were  used.  Other  life-cycle  costing  methods,  such  as 
described  in  (K'h  Life  Cycle  Costing  Instructions, 
with  anKndments"  can  also  be  applied  without  changing 
the  method  of  cost  comparison.  With  this  in  mind, 
the  example  economic  analysis  can  be  made.  Capital 
costs  will  be  considered  first. 

Capital  Costs 

Of  interest  here  is  not  the  absolute  cost  of  the 
components  of  the  solar  installation  but  the  cost 
relative  to  those  of  a  conventional  system.  Table  A I 
summarizes  the  estimated  relative  costs.  An  explana¬ 
tory  discussion  follows. 

The  values  in  the  tabic  arc  estimates  based  on 
limited  data.  The  net  cost  fur  auxiliary  heating  is 
zero,  because  it  is  sized  to  meet  the  peak  load,  the 
same  as  in  a  conventional  system. 

l.ife  Cyde  Costing  Instructions  (Department  of  the 
Army,  May  1971). 


Table  A I 

Capital  CoiU  (CC)  in  Dollan 


(.  omponent 

Cost  for 
Given  Load 

Cost/sq  ft  of 
Collector  Area 

Auxiliary  Heating  system 

0 

Cooling  System 

1548 

Controls 

3000 

Heal  Exchanger 

OiO 

Collector  l-'luid 

0.18 

Storage  T  ank 

1.00 

Plumbing 

0.75 

Pumps 

0.75 

Collectors 

5.00  or  10.00 

Labor 

1.00 

Total  CC,  =  4548  +  9.93  x  A^, 

for  less  expensive  collector  or 

CC^j  =  4548  ♦  14.93  X  A^  for  more  expensive  collector 
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h'ur  the  cooling  system  the  cost  ol  the  solar  absorption 
cooler  (including  cooling  tower)  is  compared  to  a  con¬ 
ventional,  electrically  driven  reciprocating  compretiion 
air  conditioner.  The  peak  cooling  load  (see  Table  2) 
is  4.3  X  lO’  Btu/hr  or  about  3b  tons  tor  the  Fort 
Hood  building.  Hstimated  costs  arc  $460/ton  lor  the 
gas  absorption  unit  and  $41 7/ton  for  the  compression 
unit.  This  is  a  differential  cost  of  S43/ton  or  $1548  for 
the  36  tons. 

The  cost  of  control  sensors  and  valves  is  nearly 
independent  of  system  si/.e  and  is  roughly  estimated  at 
$3000. 

Tire  rest  of  the  items  in  the  table  are  dependent  on 
the  collector  area.  Although  nut  strictly  true,  it  is 
assumed  that  the  dependency  is  linear.  Two  values  fur 
unit  collector  costs  are  used;  the  higher  value  represents 
present  low -volume  production  cost,  and  the  lower 
value  represents  a  possible  future  high-volume  produc¬ 
tion  cost. 

Fual  Coitt 

The  present-worth  cost  of  fuel  fur  the  solar  installa¬ 
tion  (PWCj)  may  be  expressed  in  terms  of  the  collector 
area  by 

PWC=PWC:L  ll-pfAJl  (KqA2) 

* 

where  p  is  the  fraction  of  the  total  energy  met  by  solar 
energy  as  expressed  in  Eq  Al  above,  and  PWC  i  a 

is  the  present-worth  cost  of  fuel  with  zero  coUector 
area. 

The  total  annual  energy  requirement  for  the  Fort 
Hood  building  (see  Table  2)  is  1 .23  X  10^  Btu.  It  is 
assumed  that  the  auxiliary  heat  is  supplied  by  gas  at  a 
current  cost  of  $l/million  Btu.  (This  was  the  June 
1975  rate  at  Fort  Hood.)  Thus,  for  Aj.»0,  the  annual 
fuel  cost  is  $1230.  Using  a  20-yr  life,  an  annual  dis¬ 
count  rate  of  10  percent,  and  an  annual  fuel  escalation 
rate  of  20  percent,  then  PWCf^^  «$68,900  and 
PWC'68.900(l-p). 

Consider  now  the  present  worth  of  the  fuel  for  the 
conventional  system  (PWC^).  Table  2  shows  that  the 
annual  cooling  load  is  7.7  X  10"  Btu.  If  the  coefficient 


of  performance  (COP)  is  3,  the  annual  electrical  con¬ 
sumption  is  7.52  X  10"  kWh.  At  $O.OI7/kWh,  the 
annual  cooling  cost  is  $1278.  Table  2  shows  that  the 
annual  heating  load  is  4  X  10^  Btu.  At  $  I /million  Btu 
fur  gas,  tlic  annual  heating  cost  is  $40.  Fur  heating  and 
cooling  tlie  annual  fuel  cost  is  $40  +  $1278  =  $1318. 
For  this  cost  and  the  same  values  fur  discount  rate, 
life,  and  escalation  rate,  it  can  be  seen  that  PWC^  = 
$76,216. 

The  net-present  worth  of  fuel  (PWC)  is  the  differ¬ 
ence  of: 

PWC  =  PWe^  PWe^  =  68,900  ( l-p|  76,216 

=  731*6  68,900  Xp. 

This  may  be  expressed  in  terms  of  collector  area  by 
using  Eq  A 1 . 

PWr=  7316  68,900X162 

X  lO"  (A  4.3  X  10  *  A 

'  c  c  ' 

=  7316  -  1 1.2  X  A  +4.8X  10  "X  A  ^ 
c  c 

The  net  life  cycle  cost  differential  (LCC)  is  simply 
the  sum  of  the  capital  cost  and  the  present  worth  of 
fuel,i.e., 

LCC,  =  CC,  +  PWC 

=  4548  +  9.93  X  A^  7316-11.2 
X  A  +4.8X  lO  "  A  ^ 

c  c 

=  -2768  -  1.27  X  A  +4.8X10"  A* 

c  c 

for  the  less  expensive  collector  ($5/sq  ft), 
and 

LCC,  =  -2768  +  3.73  X  A  =  4.8  X  10  "A  ^ 

2  C  C 

for  the  more  expensive  collector  ($10/sq  ft). 

These  two  life-cycle  costs  are  plotted  in  Figure  Al.  It 
shows  that  for  the  leu  expensive  collector,  the  solar 
installation  is  cheaper  than  the  conventional  installa¬ 
tion  up  to  a  collector  area  of  4100  sq  ft  (377  m^). 
The  savings,  although  modest,  are  greatest  for  a  col¬ 
lector  area  of  about  ISOO  sq  ft  (138  m^).  The  more 
expensive  collector  causes  the  solar  system  to  cost 
more,  although  the  added  cost  is  small  compared  to  the 
total  life-cycle  cost  of  the  system. 


33 


COLLECTOR  AREA-  10^ 


Figure  AI.  Example  of  life-cycle  costs  for  HQ  building  at  Fort  Hood,  TX. 
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APPENDIX  B; 

NBSLD-BRIEF  PROGRAM  DESCRIPTION 

The  hourly  healing  and  cooling  load  lor  each 
building  at  each  site  was  estimated  by  using  a  modiHed 
version  of  the  National  Bureau  of  Standards  l  oad 
Determining  Computer  Program  (NBSLD).  This  pro¬ 
gram  uses  the  thermal  response  factor  method  to 
estimate  conduction  and  thermal  capacitance  of  walls, 
floors,  and  roofs,  and  performs  a  detailed  hourly  heal 
balance  on  all  interior  surfaces.  The  program  pemiits 
consideration  of  variable  internal  space  temperatures 
and  provides  full  consideration  of  the  effects  of  win¬ 
dow  and  wall  shading;  internal  occupancy ,  lighting,  and 
equipment  loads;  building  use  schedules,  infiltration 
and  ventilation  air  loads;  and  building  orientation. 

Both  buildings  studied  in  this  report  were  simulated 
as  single  /.ones.  The  building's  long  axis  was  assumed  to 
be  east -west  oriented;  the  building  was  assumed  to  be 
heated  when  space  temperature  was  68“F  (25.3“C)  or 
above.  Construction  details  and  lighting  loads  were 
obtained  from  construction  drawings  of  the  buildings. 
Occupancy  and  lighting  schedules  were  estimated  from 
available  historical  data  and  information  obtained  from 
users  of  identical  existing  facilities. 

The  NBSLD  program  was  used  because  it  is  recog¬ 
nized  as  one  of  the  most  precise  load-predicting  pro¬ 
grams  available.  This  precision  was  particularly  impor¬ 
tant  to  foUow-up  solar  simulation  studies  because  of 
the  hourly  building  thermal  load's  critical  impact  on 
estimating  overall  performance  of  solar  heating  and 
cooHng  systems. 

APPENDIX  C: 

SOLAR  HEATING  AND 
COOLING  SYSTEM  SIMULATION 
PROGRAM  DESCRIPTION 

Introduction 

To  perfoim  this  study ,  it  was  nekessary  to  develop 
an  efficient  program  for  parametric  anatysis  of  solar 
heating  and  cooling  systems.  Year-long,  hour-by-hour 
system  operation  simulations  were  necessary  to  ac¬ 
count  for  daily  and  seasonal  variations  of  heating  and 
cooling  loads  and  for  the  availability  of  sunlight.  In 
order  to  save  computing  time,  the  computer  program, 
SOLSVS,  was  developed  with  a  tingle  subroutine, 
rather  than  several,  to  describe  the  entire  solar  heating 
and  cooling  system.  The  greatest  time  savings  resulted 


from  segmentation  of  the  program  and  subroutine, 
so  that  numbers  unchanged  from  the  previous  hour  or 
iteration  were  nut  recalculated.  The  data  input  method 
was  written  so  that  many  runs  could  be  made  sequen¬ 
tially,  with  only  one  or  a  few  parameters  changed  for 
each  run. 

System  Description 

Figure  Cl  is  a  diagram  of  the  simulated  solar  heating 
and  cooling  system.  The  heart  of  the  system  is  the 
storage  tank,  which  stores  energy  in  a  liquid  (probably 
water  with  suitable  corrosion  inhibitors)  in  the  form  of 
sensible  heat.  The  solar  collectors  are  separated  from 
the  storage  tank  by  a  heal  exchanger.  This  allows  use 
of  an  expansive  fluid  in  the  collectors  which  will  not 
free/e  or  boil  when  the  collector  pump  is  not  operating. 

The  heater  is  a  simple  heat  exchanger  mounted  in 
the  air  duct.  The  air  conditioner  is  a  lithium  bromide- 
water  absorption  air  conditioner.  It  was  modeled  by  a 
least  squares  fit  to  performance  curves  on  a  graph, 
showing  capacity  versus  generator  inlet  temperature 
and  outdoor  wet  bulb  temperature  for  an  ARKLA  3- 
ton  unit.*'  If  the  storage  tank  temperature  is  tix)  low 
to  operate  the  air  conditioner  or  to  provide  sufficient 
heating,  the  bypass  valve  diverts  flow  through  the  aux¬ 
iliary  heater,  which  heats  the  water  to  designated 
temperatures  for  heating  or  cooling. 

All  heat  losses,  both  from  piping  and  from  addition 
of  heat  to  the  fluids  caused  by  pump  inefficiencies, 
were  ignored.  Heating  and  cooling  loads  were  calcu¬ 
lated  separately.  The  NBSLD  program was  used  to 
create  a  data  tape  for  SOLSYS,  which  provided  hourly 
values  of  sensible  heating  or  cooling  load,  latent 
cooling  load,  dry  bulb  temperature,  wet  bulb  temper- 
aiUM,  wUkI  speed,  and  (oUd  solar  radfaUon  Incident 
on  a  horizontal  surface. 

Program  Dneription 

SOLSYS  is  cort»}jrised  of  the  main  program^  the 
system  subroutine,  and  several  small  subroutines.  The 
main  program  controls  the  data  reading,  iteration,  and 
print  processes.  Subroutine  SYSTEM  is  the  system 


’correspondence  from  Dr,  Philip  Anderson  to  Doug 
Hitlle  (CERL): Subject:  ARKLA  Cooler,  I'ebruMy  27, 1974. 

0,  Computer  Program  for  Heating  and  Cooling 
iMods  In  Building!  (National  Bureau  of  Standards,  1974). 
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Figure  Cl.  Diagram  of  simulated  heating  and  cooling  system. 


model.  Subroutines  PRNTI,  PRNT2,  PRNT3,  and 
PRNT4  print  various  outputs  called  from  SOLSYS. 
Subroutines  AC3T  and  TRNSAB,  which  are  called 
from  SYSTEM,  ate  polynomial  expressions  for  absorp¬ 
tion  8<r  conditioner  performance  and  the  transmittance- 
abs  >  rptance  product  for  the  collector  cover. 

Figure  C2  is  a  block  diagram  of  the  main  program, 
which  is  compriMd  primarily  of  an  initialization 
section  and  a  cyclic  section. 

The  initialization  section  fust  reads  the  data  cards 
which  select  the  outputs  to  be  printed  and  provide 
output  column  headings.  The  section  then  reads  a 
SCONTROL  card,  which  uses  namelist  input'*  to 
designate  the  variable  values  which  control  program 
execution.  A  check  is  then  made  on  one  of  the  control 
variables  which  continues  or  terminates  the  program. 
Next,  the  initial  values  of  the  load  and  weather  data 
and  the  initial  storage  tank  temperature  are  read  from 
an  SINITIAL  card.  Then  the  initialization  section 


FORTRAN  Extended,  Version  4,  Reference  Manual 
(Control  Data  Corporation,  1974). 


of  the  SYSTEM  subroutine  is  called,  which  reads 
parameters  describing  the  system  from  the  SPARAM 
card  and  compiites  any  values  which  are  constant 
throughout  the  fimulation. 

Next  the  program  checks  the  control  variable  which 
specifies  that  an  instantaneous  performance  calculation 
or  a  dynamic  simulation  will  be  pe 'formed.  For  a  single, 
instantaneous  calculation,  the  SYSTEM  subroutine  is 
called,  and  the  desired  outputs  are  printed.  The  pro¬ 
gram  returns  to  point  A,  where  new  control  variables, 
initial  values,  or  system  parameters  may  be  read. 

For  the  dynamic  run,  initial  values  are  set  for  all 
indexing  parameters,  and  the  data  tape  is  positioned 
to  read  the  first  day’s  data.  The  tape  reader  section 
has  been  written  for  a  specific  tape  format  and  would 
have  to  be  altered  if  a  different  format  were  used. 

Control  then  transfers  to  point  B  in  the  cyclic  sec¬ 
tion.  The  cyclic  section  also  has  a  tape  reader  to  read  a 
new  day  of  data  when  required.  This  tape-reading  sec¬ 
tion  includes  an  end-of-flle  check  to  terminate  simula¬ 
tion,  in  addition  to  an  end-of-run  check  against  one  of 
the  control  variables. 
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Figure  C2.  (cont'd). 
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The  hour  record  is  then  advanced  by  one  time 
step.  A  linear  interpolator  is  available  if  the  time  step 
is  less  than  the  time  between  data  points  on  the  tape, 
which  must  be  an  integral  multiple  of  the  time  step. 
Also,  it  should  not  be  less  than  the  time  step,  since 
some  data  wouid  then  he  ignored. 

The  day  record  is  advanced  when  hour  24  is  reached, 
and  the  solar  declination  for  that  day  is  calculated. 
Point  B  (entry  from  the  initialization  section)  is  just 
after  the  declination  calculation. 

Next,  the  iteration  section  computes  the  end  of 
time  step  tank  energy  from  the  beginning,  and  the 
SYSTtM  subroutine  calculates  the  system  performance. 
Iteration  is  terminated  when  the  relative  error  between 
two  successive  iterations  is  sufficiently  small.  There  is 
also  an  iteration  counter  which  will  temiinate  iteration 
and  the  current  simulation  by  printing  the  data  array 
for  the  day,  the  current  values  of  the  output  array, 
and  the  values  of  the  tank  energy  fur  each  iteration  in 
the  current  time  step.  These  values  can  be  helpful  for 
determining  the  reason  for  noncunvergence.  Control 
then  goes  to  point  A  for  the  next  simulation. 

Converged  iteielKi.ts  pass  to  the  integiafion  sertion, 
where  certain  outputs,  such  as  total  solar  energy  col¬ 
lected,  are  integrated  over  ume  by  means  of  a  simple, 
trapezviidat  niimerical  hitegraftm . 

Then  the  outputs  are  printed  according  to  the 
two  print  control  variables.  One  specifies  whether 
every  time  step  will  be  printed;  the  other  specifies  the 
occurrence  of  printing  once  per  day  at  midnight. 
The  daily  printer  gives  daily  integrated  values  for  all 
integrable  outputs.  The  program  then  returns  to  the 
tape  reader  section  to  begin  the  next  time  step. 

After  reaching  the  end  of  the  simulation,  a  terminal 
printout  is  made  which  includes  outputs  integrated 
over  the  entire  run.  Control  then  passes  to  point  A 
in  the  initialization  section  to  begin  a  new  simulation. 
Here,  the  value  of  namelist  input  appears;  the  nar..e- 
Ust  cards  (SCONTROL,  SINITIAL,  SPARAM)  contain 
only  variables  whose  values  differ  from  those  of  the 
previous  simulation.  All  unlisted  variables  remain 
unchanged. 

Figure  C3  is  a  diagram  of  the  SYSTEM  subroutine, 
wMeh,  h  eumpused  of  thtee  main  seetkjiis;  ittttiditea 
tkm,  control,  and  performance. 


system  parameters  and  calculates  all  values  which  will 
not  change  during  a  dynamic  simulation.  The  control 
section  is  entered  only  at  tlie  beginning  of  each  time 
step.  It  detemiines  the  mode  of  operation  of  each  sys¬ 
tem  component  and  calculates  all  values  which  are 
constant  during  the  time  step. 

The  performance  section,  which  calculates  the 
perfomiance  of  each  system  component  based  on  the 
mode  of  operation,  is  entered  during  every  iteration. 
Each  system  component  is  a  separately  written  section 
having  certain  parameters,  inputs,  and  outputs.  These 
sections  are  connected  within  the  subroutine  by 
equivalence  statements.  The  following  sections  provide 
operational  details  of  various  parts  of  the  program. 

Convergence 

The  iterative  process  is  based  on  computing  end-of- 
time-step  tank  energy  (E^),  which  is  based  on  begin- 
ning-of-time-step  energy  (E^)  and  the  change  in  energy 
caused  by  system  operation.  The  system  subroutine 
computes  rate  of  energy  change  (F  =  dE/dt)  based  on 
the  average  tank  energy  during  the  time  step  (E^  = 

t  /  2).  FXh  the  ftm  tu«l(on,  (t  (%  assumed  that 
Ej  “  ^  first  estimate  for  end-of-time-step  energy 

(E, )  is  calculated  from  E,  =  E^  +  FAt.  In  the  second 
iteTtrtvMi ,  Hm  system  idbrvwnii*  uses  t  ”  (E  +  E ,  2 
to  compute  a  new  F,  which  is  used  to  compute  a  new 
end-of-time-step  energy  (E^). 

This  iterative  process  terminates  when  the  relative 
error  between  two  successive  end-of-time-step  energies 
is  less  than  the  specified  convergence  limit  (e). 


The  significance  of  the  convergence  limit  is  more 
apparent  when  the  tank  energy  is  specified  in  terms  of 
its  temperature,  E  =  mc^T. 


T 

IVg 

c 

Maximum  value  of  T_  -  T„ 

n  n 

which  satisfies  convergence 

100  F 

.1 

10°F(-12.1‘’C) 

200  f 

.! 

100  F 

.01 

^FC-n^C) 

200  F  .01  2‘’F(-16.5‘’C) 


The  initialization  section  is  entered  only  at  the 
beginning  of  each  simulation.  It  reads  and  prints  the 
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If  the  relationship  between  tank  size  and  diange 
in  energy  is  such  that  the  tank  temperature  cannot 
change  by  the  amount  in  the  third  column,  there  will 
be  no  second  iteration,  and  the  solution  will  appear 
to  converge.  On  the  other  hand,  very  small  values  for 
tire  convergence  limit  are  unnecessary,  because  system 
inputs  are  relatively  cyclic,  and  errors  in  one  portion 
of  the  cycle  tend  to  be  canceled  by  errors  in  another 
portion.  A  limit  of  ten  iterations  was  allowed  for 
occurrence  of  convergence. 

Two  convergence  problems  were  encountered 
during  the  development  of  this  program.  The  first  was 
the  switching  between  operating  modes  in  the  SYS¬ 
TEM  subroutine,  and  the  second  was  he  model's 
inherent  instability  when  small  storage  tanks  were 
used. 

For  example,  when  the  heater  is  operational  it  will 
pump  energy  from  the  tank,  thereby  reducing  the  tank 
temperature,  as  tong  as  the  tank  temperature  is  above 
some  specified  value.  When  tank  temperature  is  below 
this  specified  value,  the  pumping  stops  and  the  follow¬ 
ing  situation  may  occur.  At  the  beginning  of  the  time 
step,  the  tank  temperature  is  just  above  the  cut-off 
value.  Energy  is  removed  from  the  tank  during  the 
first  iteration;  however,  during  the  second  iteration, 
the  average  tank  temperature  is  below  the  cut-off 
value,  so  the  pump  will  not  operate.  Ej  then  will 
equal  E^^,  since  no  energy  was  removed  from  the  tank 
during  the  second  iteration.  The  third  iteration  will 
now  be  exactly  the  same  as  the  first,  and  the  fourth 
will  be  the  same  as  the  second,  etc.  The  collector  and 
air  conditioner  have  similar  operational  cut-off  points. 

This  problem  was  handled  by  splitting  the  SYSTEM 
subroutine  into  control  and  performance  sections.  The 
control  section  determines  the  mode  of  operation  one 
time,  and  determines  whether  water  for  the  heater  is 
being  pumped  from  storage.  The  performance  Kction, 
which  is  entered  at  every  time  step,  cannot  change  the 


mode  of  operation.  This  technique  is  inaccurate  to  the 
extent  that  tlie  tank  temperature  can  change  during 
one  time  step,  and  is  thus  leu  accurate  fox  smaller 
tanks;  however,  the  cyclic  nature  of  the  dynamic  prob¬ 
lem  causes  error  cancellation  over  a  lung  time  period. 

When  small  storage  tanks  were  used,  several  cases 
of  nonconvergence  were  observed  where  successive 
iterations  oscillated  around  an  apparent  convergence 
point  (Sec  Figure  C4).  The  convergence  was  too  slow 
for  the  convergence  criterion  to  be  met  in  ten  iterations. 


By  defining 


Al  -  E ,  -  E 

I  <1 


A2  =  E,  E, 


A3  =  E3  E, 


It  was  observed  that  A1/A2  s  A2/A3  a  A3/A4  . . . 

...  *  R,  it 


Assuming  A I /A2  «  A2/A3  =  A3/A4 
can  be  shown  that 

Ea-E| 

E  =  E,  - - 

*  *  1  -R 


(EqCl) 


This  method  succeufully  extended  the  lower  limit  of 
tank  size  to  about  one-half  of  its  original  value.  Further 
tank  size  reduction  led  to  divergence,  so  the  assump¬ 
tion  of  constant  ratios  was  no  longer  valid;  however, 
convergence  could  still  be  achieved  by  reducing  the 
time  step  size. 

It  was  considered  desirable  to  use  the  longest  time 
step  possible  to  minimize  iterations  for  a  l-yr  run. 
Nonconvergence  with  the  corrector  occurred  when 
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Figure  C4.  Energy  vs.  number  of  iterations  for  slowly  converging  case. 
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(i;qC5) 


lank  d/e  was  so  small  that  a  20-degree  change  in  tank 
temperature  occurred  at  maximum  load.  Such  an 
extreme  change  contradicts  the  quasisteady  assumption 
used  to  derive  the  system  performance  equations. 
Therefore,  it  is  probably  not  desirable  to  pursue 
additional  means  of  achieving  convergence  without 
reducing  the  time  step. 

Solar  Radiation  Equations 

Solar  radiation  is  input  as  the  intensity  (Btu/hr/ 
sq  ft)  on  a  horizontal  surface.  This  figure  must  be 
converted  into  components  of  beam,  diffuse,  and 
reflected  radiation,  so  that  collector  performance  can 
be  calculated.  This  was  done  by  following  the  tech¬ 
nique  presented  by  ASHRAb  and  used  in  NBSLD.'^ 

ASIIRAb  gives  the  intensity  of  beam  radiation  as 

1^  =  Ae<-®/“"'’>  (BqC2) 

where  A  =  apparent  solar  irradiation  at  air  mass  0 
B  <■  atmospheric  extinction  coefficient 
0  »  solar  altitude 

The  diffuse  radiation  on  a  horizontal  surface  is  given  by 
lj=CI,  (EqC3) 

where  C  >=  diffuse  radiation  factor 

NBSLD  modifies  these  two  equations  to 

I'  -  DAe<-«/*'"^) 
and 

I'  »  CI'/D* 

where  D  =  atmospheric  clearness  number'  ^ 

From  these  two  equations,  the  intensity  of  radiation 
on  a  horizontal  surface  is 

‘h  “  +  (EqC4) 


^^ASHRAE  Handbook  of  Fundamentalt  (American Society 
of  Healing,  Refrigerating,  and  Air  Conditioning  F.ngiiwcri, 
Inc.,  1972),  p  394. 

'^J.  A.  Duffle  and  W.  A.  Beckman,  Sbfer  Energy  Thermal 
Procetut  (John  Wiley  and  Sons,  1974),  p  15. 


Therefore,!^  *  l^/(sin/}  +  (7D*). 

The  angle  of  beam  radiation  incidence  on  a  tilted 
surface  is  given  by 

cos  0^  =  sin  h  sin  0  cos  s  sin  A  cos  it>  sin  s  cos  7 
-f  cos  A  sin  0  cos  s  cos  co 
4  cos  A  sin  0  sin  s  cos  y  cos  co 
-r  cos  A  sin  s  sin  7  sin  a;  (l-.qCo) 

where 

0  =  latitude  of  site  (north  positive) 

A  =  solar  declination 
s  *  slope  of  the  surface 
7  =  azimuth  of  the  surface  (south  =  0,  east  =+) 
CO  -  hour  angle 

The  angle  of  incidence  on  u  horizontal  surface  is  given 
by 

cos  0^  =  sin  A  sin  0  +  cos  6  cos  0  cos  w  (tq  C7) 
Note  that  cos  ^  =  sin  0. 

The  view  factor  between  the  collector  and  the  sky  is 
=  (I  ■•■cos  s)/2.  The  view  factor  between  the 
collector  and  the  ground  if  F^,^  =  ( 1  cos  s)/2 . 

The  total  radiation  incident  on  the  collector  is 

1|  =  beam  radiation  ■•■diffuse  radiation  ■•■  reflected 
radiation 

-l^cosOc  +  l'  F„  +  (lhCosfl^  +  l')p,F^^ 

(EqCS) 

where  =  reflectivity  of  the  ground. 

To  avoid  recalculation  of  the  sine  and  cosine  terms, 
the  radiation  calculations  appear  in  several  parts  of 
the  program.  Since  latitude,  slope,  and  azimuth  will 
not  change  during  a  simulation,  they  are  calculated 
once  in  the  initialization  section  of  SYSTEM.  Dedina¬ 
tion  and  the  diffuse  radiation  factor  are  computed 
once  daily  in  the  main  program.  Angle  and  radiation 
intensity  are  computed  once  hourly  for  each  time 
step. 
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Colltctor  Equatiom 

Figure  CS  illuttratet  the  lyitem  configuration  and 
variables  uied  in  the  collector  equationi. 

Parameters; 

A  collector  area 

(ra)  iransmittance-absorptance  product 

^pc  ipcciflc  heat  of  collector  fluid 

mau  flow  rate  of  collector  fluid 

F'  collector  geometry  efficiency  factor'^ 

edge  and  back  heat  Ion  coefficient 

N  number  of  glass  covers 

collector  plate  emiuivity 
F  heat  exchanger  effectivenesa 

Cp,  specific  heat  of  tank  fluid 

tfi|  mass  flow  rate  of  tank  fluid 

Inputs: 

S  solar  radiation  on  the  collector 


ambient  dry  bulb  temperature 
w  wind  speed 
T. 


I  temperature  of  fluid  from  tank 
Computed  values: 


til  f  ' 

t  pc 


C  =rh  C 

t  c  pc 


c 

c 

C  ,  =  min  of  — 
min  P 


U|  A 

il 

fl  e  ^ 


F.C 


P‘ 


min 


The  rate  of  energy  collection,  q^,  can  be  expressed 
by  the  following  formula: '  ’ 

q^*AF,|S(Ta)-UL(T^ -TJl  (EqC9) 

Also,  q^  «  -  T^)  or  T,  »  +  qJC^ 


(EqCIO) 


The  heat  exchanger  equation  is 


Ta  = 


FCmJT.-T,) 


+  T,  (EqCIl) 


'^R.  W.  Htiss,  "The  Derivition  of  Several  Plate  Efficiency 

l  actori  UKful  in  the  Deii|n  of  l  lat'Plalc  Solar  Heal  CoUec-  A.  Duffle  and  W.  A.  Beckman,  So/ar  Energy  Thermal 

torj,"  Solar  Energy,  Vol  3,  No.  4  (1959),  pp  55-64.  Procenet  (John  Wiley  and  Sons,  1974),  p  241 . 


Figute  CS.  CoUector-ttorage  tank  system  schematic. 
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There  are  three  unknowns  (q^,  T^,  Tj )  in  these  three 
equations.  Substituting  Kqs  CIO  and  Cl  I  into  Kq  C‘) 


{f/ISIto)  U,  ( 

t 

Tjl).  - ; 


T,) 


«.in<T,  T,) 


F/siro)  f;au,_f;c„,j, 


N  I 

C  (T  T)  ”  h 

H - 

T  N+l 


(T  +T  KT  +T  ) 

■ - +  U, 

I  2N  +  f  I  '’*■ 

- + - —  N 

e^  +  .OSNd  fp) 


(KqC  1.3) 


|KJ  hr  '  Ml  ^  K  'I 


F  All,  FC  ,  T,  F' AU,  T, 

r  1,  min  2  r  L  2 


FAU,  T  f;c  .  t, 

r  L  V  min  I 

+  - + - 

c  c 

c  c 


f:c  t, 

m  in  2 


- +  T- 

C.  * 


Tp  =  mean  collector  plate  temperature 
=  (Tj+T^)/2 


h  =  20..S2+  13.68  w 

W 


f  =  (I  .04h^  +  .OOOSh^KI  +  .0‘)IN) 
C  =  365.9(1  -  .00883s +.0001 298s^) 


oS(Ta) 

Tj  = -—  Q/3T, +a^Tj  oTj+aT^ 
^1. 


+  0T,  -/3Tj  +  1^ 


aS(ra) 


+  (0-o0)T,  +oT 


s  =  collector  slope  (degrees) 

Since  U,  is  a  function  of  T,  and  T  .  its  value  is 

L  2  4 

most  conveniently  found  iteratively;  however,  since 
it  is  a  weak  function  of  and  T^,  it  is  evaluated  only 
on  the  First  iteration  in  each  timestep. 

The  transmittance-absorptance  product  is  given  by 


a  +  0  -  a/J 


The  other  unknowns  can  now  be  found. 


(EqCI2) 


1  -  0  -  a)p . 


(EqCI4) 


%  =  EC„JT,-T,) 

T.=T,-q„/C^ 

T3=T,+q„/C, 


The  overall  heat  loss  coefficient,  U^,  is  a  function 
of  collector  construction  and  operating  conditions. 
The  following  expression  was  developed  by  S.  A. 
Klein' ^  from  experimental  studies. 


a  =  collector  plate  absorptance 

1  -p 

T  =  cover  transmittance  = - X  e  *' 

1  +(2N-I)p 

^  r  due  to  reflection  *  r  due  to  absorption 
,  sin’(0j-fl,)  tan*(0j-^,) 

p  *  —  [ - + - ) 

2  sin*(0j+fl,)  tan*(flj+0,) 

where  0 1  and  are  as  shown  in  Figure  C6 


"*S.  A.  Kkin,  “Calculation  of  Flat-PUte  Collector  Loss 
Coefndcnti,”  paper  presented  at  U.  S.  Section  meeting, 
International  Solar  Energy  Society,  Fort  Collins,  CO  (1974). 


''J.  A.  Duffle  and  W.  A.  Beckman,  5»/0r  A'neryv  Thermal 
Provettes  (John  Wiley  and  Sons,  1974),  pill. 


tcmpcruturc  from  auxiliary  heater 


K  =  extinction  coefficient 

L  =  distance  light  travels  through  the  cover. 

To  save  computer  time,  a  least  squares  fit  was 
performed  for  I.  2,  and  3  cover  systems  and  the  re¬ 
sults  expressed  as  a  sixth-order  polynomial.  The  poly¬ 
nomial  is  accurate  to  within  I  percent  fur  low  iron 
glass. 

The  collector  has  two  modes  of  operation:  both 
pumps  either  may  be  pumping  or  not  pumping.  This 
model  assumes  pumping  whenever  the  energy  collected 
is  positive.  Actual  operation  would  involve  temperature 
sensors  which  would  start  the  pumps  whe.:  Tj  <  Tj. 


Figure  C6.  Incident  and  refracted  angles. 


Haatar  Equationi 

The  healer  shown  in  Figure  C7  is  operated  by  either 
stored  or  auxiliary  heal.  Actual  control  of  this  heater 
would  involve  a  dual  deadband  with  three  tempera¬ 
tures;  T^>T  >T^.  When  the  room  air  temperature,  T^, 
falls  below  the  pump  starts;  if  T^ rises  above  T^, 
the  pump  stops.  But  if  T^  faDs  below  T^,  the  valve 
will  divert  the  flow  through  the  auxiliary  heater,  whidi 
will  operate  until  T  ,>Tj,.  The  auxiliary  heater  operates 
by  heating  the  water  to  a  specified  temperature,  T,. 
If  the  storage  temperature  it  to  low  (T,<T|jj||j)  that 
it  is  inefficient  to  pump  from  the  storage  tank,  the 
valve  will  always  divert  through  the  auxiliary  heater. 

Parameters: 

m^  mass  flow  rate  produced  by  pump 
fluid  specific  heat 

ECm,„  heat  exchanger  effectiveneu  multiplied  by 
minimum  heat  rate  (heat  rate  of  the  air) 

^min  minimum  storage  tank  temperature 


T, 

T^  mom  air  tcmpcraturc(treatcdasaconstant) 

Inputs: 

T|  lank  temperature 

‘lio.d  I'ca'inKl'wd 

Outputs: 

heating  load  met  from  storage 

q^  heating  load  met  from  auxiliary 

Tj  temperature  of  fluid  returning  to  storage 

lank 

th  average  mass  flow  rate  to  and  from  storage 

The  equations  describing  this  system  are  based  on 
two  maximum  heating  rates  and  the  heat  exchanger 
equation.  The  maximum  rate  of  transferring  heat  from 
the  storage  tank  to  the  room  air  is 

(q  )  =  KC  .  (T,  -  T  )  «  m  c  (T,  -  T, ) 

(F.qC15) 

The  maximum  rate  of  transferring  heat  from  the  auxil¬ 
iary  heater  is 

=  f^C'^,„(T3-T.)  =  m.c;(T3-T,) 

(EqC16) 

This  heater  has  four  modes  of  operation.  First  T| 
may  be  less  than  T  ,  ;  in  this  case, m  =  0,  q  =0,and 
4  =  <li  Note  that  (q  )  must  be  greater  than 
the  maximum  probable  heating  load.  This,  in  combina¬ 
tion  with  the  auxiliary  heater  temperature,  T^,  deter¬ 
mines  the  size  of  the  heat  exchanger. 

In  the  second  mode,  the  entire  heating  load  can  be 
met  by  pumping  from  storage;  that  is,  when  (q  ) 

‘fio.d’  =  ‘lioid  ‘I,  “  «  ''-hou* 

perM,  the  entire  heating  load  can  be  met  by  pump¬ 
ing  for  only  a  fraction  (Y)  of  the  hour,  where  Y  = 
‘l|iiad/^‘'ii^mii»‘  average  flow  rate  for  the  stor- 
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Figure  C7.  Auxiliary  heater  circuit . 


age  tank  is  rii  =  Y  iti^,  and  during  pumping,  Tj  =  T, 
(q  )  /rii  c  . 

In  tlie  third  operating  mode.  (q,)m,„<qi„,d-  Suffi¬ 
cient  heat  must  then  be  added  from  the  auidtiary  to 
meet  the  load.  To  make  maximum  use  of  stored  en¬ 
ergy,  heat  will  be  transferred  from  storage  whenever 
the  auxiliary  is  nut  operating.  If  Y  is  the  fraction  of 
each  hour  that  heat  is  being  transferred  from  storage, 
then, 


(EqCI7» 

^  ^loid 

K>n,.x-K)mlx 


The  values  of  the  outputs  are  th  =  Y  th^,  T,  *  Tj  - 
(q  )  /rh  c  ,  q  =  Y  (q  )  ,  and  q  =  ( I  -V^q  ) 


Finally,  in  the  fourth  mode  there  may  be  no  heating 
load.  Then  m  *■  *  0. 

In  this  analysis,  the  thermal  capacity  of  the  heater 
circuit  fluid  hu  been  neglected.  This  should  be  a 
reasonable  anumption  if  switching  between  auxiliary 
heating  and  storage  heating  occurs  only  a  few  times 
per  hour,  and  there  is  relatively  little  piping  between 
the  auxiliary  and  the  heat  exchanger. 


Absorption  Air  Conditioner  Equations 


low-temperature  heat  as  its  primary  energy  source. 
Solar  air  conditioning  permits  greater  year-round  use 
of  solar  collector  and  storage  components.  Absorption 
air  conditioners  currently  in  production  have  not  been 
optimised  to  operate  at  temperatures  produced  by 
solar  collectors.  The  following  model  is  based  on  a 
3-ton  (36,000  Btu/hr)  LiBr-water  unit  designed  to 
operate  using  water  at  2I0°F  (93.2°C),  but  capable  of 
using  it  at  temperatures  as  low  as  170°F  (76.6“C) 
before  crystallization  occurs. 

The  air  conditioning  system  shown  in  Figure  C8 
uses  the  same  type  of  circuit  as  the  heater,  so  its  oper¬ 
ating  equations  are  similar. 

Parameters: 


>*0 

mass  flow  rate  produced  by  pump 

fluid  specific  heat 

CAP 

air  conditioner  capacity 

COP 

coefficient  of  performance 

T 

'min 

minimum  generator  inlet  temperature 

(T  ) 

maximum  ambient  temperature  for  no 
cooling 

Inputs: 

T, 


An  absorption  air  conditioner  is  especially  attrac¬ 
tive  for  tolar  energy  cooling  because  it  uses  relatively 
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storage  tank  temperature 


cooling  load 


and  the  maximum  cooling  rate  using  auxiliary  energy  is 


ambient  dry  bulb  temperature 
ambient  wet  bulb  temperature 

The  energy  used  to  obtain  a  given  amount  of  cooling  is 
Outputs;  Q  , 

0,  energy  fromstorage  used  for  air  conditioning 

There  are  four  operating  modes.  In  the  first  mode, 

0,  energy  from  auxiliary  used  for  air  condi-  Then  m  “0, 0^  *0  , and  0^  =  q|„^j/COP. 

tioning 

Tj  temperature  of  fluid  returning  to  storage 
tank 

rh  average  mess  flow  rate  to  and  from  storage. 

The  air  conditioner  absorber  and  condenser  are 
cooled  by  a  cooling  tower.  The  air  conditioner  per¬ 
formance  is  expressed  as  a  function  of  rated  capacity, 
based  on  the  generator  inlet  and  outside  wet  bulb 
temperatures;  therefore,  the  maximum  cooling  avail¬ 
able  by  using  stored  energy  is 

K)n,.x*CAP.f(T,.TJ 


TO 

COOLING 

TOWER 


FROM 

COOLING 

TOWER 


In  the  second  operating  mode,  the  entire  cooling 
load  can  be  met  from  storage.  0,  ®  and  0 

=  0.  The  fraction  of  each  hour  requiring  pumping  is 

mc  =  Ym  andT2  =  TI 

m^c^COP 

In  the  third  mode,  (q,)mgx^i,„d-  Letting  Y  be 
the  fraction  of  each  hour  that  energy  is  removed  from 
storage, 

Y  *  - 


Figure  C8.  Solar-powered  air-conditioning  system  schematic. 


The  values  of  the  outputs  are 


fii  =  Yrii  .  T  =  T, 

c  2  I 


th  c  COP  ■ 

c  p 


0  =  Y(q  )  /COP.  and  0  = 

s  mu*  u 


(1  YKq.)„ 


Kinaily,  there  may  be  no  cooling  load.  It  is  also 
assumed  that  there  is  no  energy  demand  when  there  is 
a  cooling  load,  but  that  the  outside  temperature  Is 
below  some  specified  value,  e.g.,  60°F  ( 1 5.5"C|.  Then, 
cither  the  cmiling  system  will  include  means  of  intro¬ 
ducing  outside  air,  or  the  building  occupai  '.swill  open 
windows. 


Storage  Tank  Equatktns 

Energy  is  stored  as  sensible  heat  in  fluid  in  the 
storage  tank  and  is  lost  from  the  storage  tanks  to  its 
surroundings.  The  loss  coefficient,  U,  is  a  parameter 
which  must  be  calculated  for  the  type  and  amount 
of  tank  insulation.  The  model  assumes  a  constant 
value  for  the  temperature  of  the  surrounding  environ¬ 
ment.  (While  this  would  be  a  good  assumption  for  a 
tank  inside  the  house,  it  is  probably  not  a  good 
assumption  for  the  tank  buried  outside  the  house;  no 
better  model  was  readily  available,  however.)  Only  one 
card  would  have  to  be  changed  for  the  energy  Ion 
from  the  tank  to  be  subtracted  from  the  heating  load 
or  added  to  the  cooling  load. 

It  is  possible  that  the  temperature  at  which  energy 
is  collected  will  exceed  the  boiling  temperature  of  the 
storage  fluid.  The  model  anumes  that  boiling  may 
occur,  with  energy  being  dissipated  through  a  preuure 
relief  valve.  The  boiling  temperature  is  a  system  para¬ 
meter. 


The  energy  stored  in  the  tank  is  given  by  E  =  mc^T, 
with  T  =  0  defining  the  zero  energy  level.  The  rate  of 
tank  energy  change  is  the  sum  of  the  energy  rates 
produced  or  consumed  by  the  other  components. 

Rate  of  energy  gain  from  collector 

q,  =  n\.Cp(T,  T,) 

Kate  of  energy  loss  to  heater 

"'hS-O,  ■rh> 

Rate  of  energy  loss  to  air  conditioner 
Rate  of  energy  loss  to  environment 

q,=  i'ait,  t.) 

where  rii^, and  rh^  arc  the  average  flow  rates  to  the 
components,  and  T^,  T^,and  T^  are  the  return  temper¬ 
atures.  T|  is  the  average  tank  temperature,  and  A  is 
the  tank  surface  area,  yielding 

dE 

M,  I,  M, 

The  safety  relief  compares  the  tank’s  total  energy 
to  the  maximum  energy  permitted  by  the  boiling 
point  limit. 

E  =  E ,  +  F  A  t 

Energy  lost  by  boiling  in  the  tank  is 
q  =L  E  ifE>E 

^Dver  mu*  max 


if  E< 


E 

max 


This  gives  a  revised  rate  of  change  of  tank  energy. 


^iver 

F  =  F- -  (EqCI9) 

At 
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